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I.  INTRODUCTION 


During  the  period  1963-1966,  Dr.  P.  C.  Plen  of  the  Naval  Ship  Research 

and  Development  Center  conducted  an  Investigation  Into  the  application  of 

12  3 

existing  wavemaking  resistance  theory  to  the  design  of  practical  ship 

hull  forms.  During  the  three  years  of  the  study.  Dr.  Plen  developed 

several  computer  programs  which  were  later  used  In  designing  the  hulls  of 

various  ship  models.  The  procedures  developed  and  the  resistances  obtained 

4  5  6 

were  documented  by  Plen.  Although  the  models  built  from  these  designs 
had  low  resistances,  the  agreement  between  the  values  of  the  experimental 
and  theoretical  resistances  was  unsatisfactory. 

With  the  recent  surge  of  Interest  In  unconventional  hull  forms,  such 
as  catamarans,  small  waterplane  area  twin-hull  (SWATH)  ships  and  small 
waterplane  area  single-hull  (SWASH)  ships  for  which  little  or  no  experi¬ 
mental  Information  Is  available,  the  theoretical  approach  to  ship  design 
has  assumed  a  new  Importance.  Accordingly,  the  programs  developed  earlier 
by  Dr.  Plen  have  been  redesigned  Into  a  new  program  called  "The  Plen  Wave¬ 
making  Resistance  Computation  Program."  This  program  computes  and  prints 
out  the  wavemaking  resistance  coefficients  (Cw)  for  various  ship  speeds. 

It  can  also  compute  the  frictional  resistance  coefficients  (Cf)  and  there¬ 
fore  the  effective  horsepower  (EHP)  for  various  ship  speeds.  Using  this 
program,  the  designer  is  able  to  predict  the  resistance  of  a  SWATH  or 
SWASH  ship  of  a  particular  design  without  undergoing  the  expense  of  build¬ 
ing  and  testing  a  scale  model. 

The  Plen  program  has  been  used  successfully  in  the  past  few  years  for 
the  design  of  several  SWATH  ships,  models  of  which  have  actually  been 
built.  Experimental  results  have  shown  close  agreement  with  the  resis- 

O 

tances  predicted  by  the  Plen  program.  (See  also  Appendix  B.) 

The  Plen  program  consists  of  about  60  individual  subroutines  assembled 
Into  three  main  overlays  which  are  available  either  on  disk  or  on  magnetic 
tape.  These  routines  are  coded  In  CDC  Fortran  Extended  for  use  on  the  CDC 
6000-series  computers.  The  program  requires  a  minimum  of  35K  octal  core 

T - 

A  complete  list  of  references  is  given  on  page  149. 
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for  execution.  An  Interactive  graphics  capability  operable  on  the  CDC 
6000/1700/274  computer  configuration  Is  being  developed  as  a  separate 
program. 

Documentation  of  the  Plen  program  has  been  prepared  In  two  parts. 

The  first  describes  the  theoretical  basis  and  programming  technique,  and 
is  Intended  primarily  for  those  who  wish  to  review  the  theoretical  equa¬ 
tions  or  understand  the  coding  of  the  program.  The  second,  the  User's 

q 

Manual,  Is  Intended  for  those  who  are  responsible  for  designing  a  ship 
hull  with  low  wavemaking  resistance  or  for  computing  the  wavemaking 
resistance  of  a  ship  hull  for  which  the  singularity  distribution  is  known. 
Programming  knowledge  Is  not  essential  to  understand  the  Users  Manual. 


II.  THEORETICAL  BASIS 


Normal  pressures  vary  along  the  length  of  a  deeply  submerged  body 
traveling  horizontally  at  a  steady  speed.  If  the  fluid  surrounding  the 
body  is  nonviscous,  the  net  force  of  these  various  pressures  is  zero. 

However,  when  a  body  travels  on  or  near  the  surface,  the  variations  in 
pressure  along  the  length  of  the  body  cause  waves  to  form  which  themselves 
affect  the  distribution  of  the  pressures  along  the  hull  so  that  a  net 
fore-and-aft  force,  called  wavemaking  resistance,  is  created.  A  great 
deal  of  research  has  been  devoted  to  the  development  of  theoretical  methods 
for  calculating  wavemaking  resistance  and  to  their  experimental  verifica¬ 
tion.  One  such  method  involves  determing  the  flow  around  the  hull,  and 
hence  the  normal  pressure  distribution,  and  then  integrating  the  fore-and- 
aft  components  of  these  pressures  over  the  hull  surface.  This  method  was 
developed  in  1898  by  Michel  1  for  a  thin  ship  moving  over  the  surface  of  a 
nonviscous  fluid  J  Another  method  is  that  used  by  the  Men  program  which 

adapts  Havelock's  mathematical  representation  of  the  wave  system  generated 

2  3 

by  the  ship  at  a  great  distance  astern  to  compute  the  wavemaking  resis¬ 

tance  as  determined  by  the  flow  of  energy  necessary  to  maintain  the  wave 
system. 

A  ship  hull  moving  through  water  can  be  hydrodynamically  represented 
by  singularity  distributions.  The  basic  element  of  a  singularity  distri¬ 
bution  is  a  simple  point  source.  From  It, a  point  doublet  can  be  derived. 

Since  a  line  source  and  doublet  and  a  surface  source  and  doublet  can  all 
be  expressed  in  terms  of  a  point  source,  only  those  equations  of  amplitude 
function  and  induced  velocity  components  of  a  point  source  are  given. 

WAVE  AMPLITUDE  FUNCTIONS  OF  A  POINT  SOURCE 

The  gravity  waves  created  by  a  point  source  with  a  steady  motion 
under  a  free  surface  have  two  parts: 

c(x,y)  =  cw(x,y)  +  ^(x.y)*  (1) 

*The  coordinate  systems  are  assumed  to  move  with  the  ship,  with  the  z-axis 
vertically  upward,  z  =  0  on  still-water  surface,  and  the  x-axis  opposite 
to  the  ship's  course. 
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where  cw(x,y)  Is  the  free-wave  pattern 
tl(x,y)  Is  the  local  disturbance 

The  free-wave  pattern  Is  the  cause  of  wave  resistance;  the  local  distur¬ 
bance  has  nothing  to  do  with  the  transfer  of  the  wave  energy. 

Sir  Thomas  H.  Havelock's  papers  treated  a  free-wave  pattern  as  a 
resultant  of  the  "elementary  waves,"  and  obtained  the  following  extremely 
simple  and  elegant  relationship  between  the  elementary  waves  and  the  wave¬ 
making  resistance.  For  a  moving  point  source  located  at  a  point  U,  n,  t)» 
the  free  surface  wave  elevation  can  be  expressed  in  the  x-y-z  coordinate 
system  as 

tw(x,y)  *  Y~  U°  fit  sec3e  *  e"k°  ^  ’  sec  9  •  cos  R  •  de  (2) 

"2 

o 

where  R  =  k  sec  e  •  [(x  -  5)  •  cos  0  +  (y  -  n)  •  sin  0] 

k0  -  g/v2 

M  *  Source  strength 
U  =  Speed  of  advance 

0  =  Angle  between  wave  direction  and  x-axis 

a  =  tan-1  (-x/y) 

Introducing 

RQ  =  X  •  cos  0  +  y  •  sin  e  (3) 

S  =  £  •  cos  0  +  n  •  sin  8  (4) 

we  can  write 

R  =  kQ  •  sec20  •  (Rq  -  S)  (5) 
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(6) 


f 


and 


cos  R  =  cos  (k  •  secze  •  S)  •  cos  (k  •  secZ9  •  R  ) 

0  0  o 

p  p 

+  sin  (kQ  •  sec  e  •  S)  •  sin  (kQ  •  sec  e  •  RQ) 
From  equations  (2)  and  (6)  we  have 

cw(x,y)  =  L  Ac(e,  c.  c)  *  cos  (KQ  •  sec2e  •  RQ)  •  de 


+  /,  A$(e,  5.  s)  •  sin  (KQ  •  secZe  •  RQ)  •  de  (7) 


where  Ac(e,  5,  c)  = 


M  •  K  -  k  .  c  •  sec  e  9 

tt~-  ~U°~  '  se<-  6  *  e  °  ‘  cos(kQ  *  sec0  •  $) 


V*.  .  u 


M  •  k  -  k  •  c  •  sec  e  9 

- - -  •  secJe  •  e  0  •  sin(kQ  •  sec^e  •  S) 


Since  RQ  Is  the  distance  between  the  point  (x,y)  and  the  origin,  the 
free-wave  system  r  (x,y)  as  shown  by  Equation  (7)  consists  of  two  com- 

W 

ponents:  namely,  a  cosine-wave  system  with  an  amplitude  function  Ac(e,  5,  5), 
and  a  sine-wave  system  with  an  amplitude  function  A  (e,  5,  c). 

The  computation  of  the  wavemaking  resistance  for  a  surface  or  a 
line  singularity  distribution  Is  based  on  the  results  of  a  point  singu¬ 
larity  moving  underneath  a  free  surface.  For  a  moving  point  source 
located  at  a  point  U,  n,  c),  a  three-dimensional  surface-wave  system  is 
created. 

Havelock  obtained  the  following  general  expression  for  wavemaking 
resistance  In  regard  to  a  point  source: 

/■»/ 2 

Rw  s  ir*p*u2  Jo  [Ac(e,s,c)]2  +  [A S(e,s,c)]2  •  cos3e  •  de  (8) 
where  p  Is  the  water  density. 
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VELOCITY  COMPONENTS  OF  A  POINT  SOURCE 


When  the  singularity  distribution  representing  a  hull  form  Is  known, 
the  hull  geometry  can  be  obtained  by  tracing  a  number  of  streamlines. 

These  streamlines  define  a  closed  stream  surface  to  approximate  the  hull 
geometry.  Each  streamline  Is  determined  from  a  set  of  first-order  dif¬ 
ferential  equations 

dy  _  v 

dx  '  U  +  u 

dz  _  w 

Hx  '  U  +  u 

where  U  Is  the  uniform  stream  velocity,  and  u,  v,  and  w  are  the  induced 
velocity  components  resulting  from  a  given  singularity  distribution.  Since 
the  induced  velocity  components  of  a  given  singularity  distribution  can  be 
derived  from  point  sources  alone,  only  the  expressions  needed  for  computing 
the  Induced  velocity  components  for  a  point  source  will  be  given  here. 
Theoretically,  the  Induced  velocity  components  of  a  point  source  should 
take  Into  consideration  the  free  surface  effects.  However,  at  the  present 
time  It  Is  only  practical  to  trace  streamlines  under  the  assumption  that 
the  free  surface  Is  a  rigid  wall.  This  Is  equivalent  to  assuming  a  double 
model  In  an  infinite  fluid  without  a  free  surface.  Under  the  rigid  wall 
condition,  another  point  source  Is  Introduced  as  the  mirror  image  for  each 
point  source  under  a  free  surface.  Thus,  the  velocity  potential  of  a 
point  source  under  a  free  surface  can  be  derived  as  follows: 

,  =  _JL_ 

4n  •  r 

where  M  =  point  source  strength 


III.  APPL I CATION  OF  THE  THEORY 


Michel 1 ' s  thln-shlp  wavemaking  resistance  theory  rests  on  two  basic 
assumptions  that  limit  Its  usefulness  for  the  design  of  practical  hull 
forms:  The  first  Is  that  the  ship  is  thin  and  the  free-surface  wave 
slope  Is  small;  the  second  Is  that  no  viscosity  effects  exist.  Although 
no  satisfactory  method  for  dealing  with  viscosity  effects  has  yet  been 
devised,  ways  have  been  found  for  dealing  with  the  thln-shlp  assumption. 

For  a  thin  ship,  the  flow  over  the  ship's  surface  Is  essentially 
two-dimensional,  so  that  the  hull  form  may  be  represented  by  a  source 
distribution  on  a  central  plane  and  the  density  of  the  source  distri¬ 
bution  may  be  assumed  to  be  proportional  to  the  waterline  slope  at  various 
depths.  In  reality,  however,  a  practical  hull  form  must  have  a  beam- 
draft  ratio  greater  than  2,  and  the  flow  around  it  will  be  highly  three- 
dimensional. 

Inul^  used  an  Inverse  method  In  which  known  singularity  distribu¬ 
tions  were  used  as  the  starting  point  and  the  corresponding  hull  form 
geometry  was  obtained  by  tracing  streamlines  of  the  singularity  distribu¬ 
tions.  These  same  singularity  distributions  were  also  used  to  compute 
the  wavemaking  resistance.  He  obtained  close  agreement  between  experi¬ 
mental  results  and  theoretical  predictions  by  building  his  experimental 
model  to  conform  with  the  steam  surface.  The  Plen  program  uses  Inul's 
Inverse  method  technique. 

However,  Inul's  approach,  too,  has  shortcomings:  The  hull  form  Is 
nearly  as  Impractical  as  Michel  1 's— due  to  the  small  beam  draft— and  the 
streamline  tracing  must  be  performed  under  zero  Froude  number  conditions. 
Although  the  Plen  program  must  also  operate  under  the  zero  Froude  number 
conditions,  since  tracing  streamlines  at  nonzero  Froude  number  conditions 
is  extremely  difficult,  the  problem  of  the  small  beam-draft  ratio  has 
been  solved.  The  Plen  program  uses  a  subsurface  (n-surface)  rather  than 
a  central  plane  as  the  location  of  the  singularity  distribution.  This 
technique  Is  explained  In  the  following  sections. 


:;v^'  P"  ' 


n-SURFACE  AND  COORDINATE  SYSTEMS 


In  the  past,  the  singularity  distribution  has  customarily  been 
located  on  the  centerline  plane  of  a  ship,  a  fact  which  has  resulted  in 
the  beam-draft  ratio,  B/T,  always  being  less  than  2.0.  To  obtain  more 
practical  B/T  values— values  greater  than  2.0— the  Plen  program  defines 
the  singularity  distribution  on  an  n-surface  as  follows: 

n  =  +  f(c.  c)  (11) 


Figure  1  -  Coordinate  Systems  and  n-Surface  for 
Location  of  Singularity  Distributions 

where  5,  n.  and  c  are  the  coordinates  of  a  rectangular  coordinate  system 
with  the  origin  at  the  midship  on  the  undisturbed  free  surface. 

Let's  define 

n  =  Bn  •  (1  -  cn)  (12) 

where  n  Is  an  even  Integer  number  and  Equation  (12)  represents  a  vertical 
strut-like  surface,  Its  beam  controlled  by  B^  and  Its  waterline  controlled 
by  n. 

To  facilitate  numerical  Integration,  a  waterline  of  an  n-surface  is 
further  simplified  by  being  represented  as  a  number  of  first-degree 
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The  n-surface  Is  then  represented  by  K  piecewise  straight-line  segments  as 


where  k  =  0,  1,  2  •••  K 

In  the  computer  program,  K  Is  arbitrarily  set  less  than  5,  1  e., 

1  s  k  s  4. 

By  varying  the  parameters  of  the  n-surface  and  the  singularity  distri¬ 
bution,  a  large  family  of  different  hull  forms  can  be  generated.  By 
choosing  the  appropriate  beam  B^  and  depth  of  the  n-surface,  the 
required  B/T  value  of  a  set  of  ship  hull  forms  can  be  obtained. 
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SINGULARITY  DISTRIBUTIONS 


The  singularity  density  distribution,  defined  as  the  singularity 

M 

strength  per  unit  area  per  unit  velocity  of  a  moving  ship,  jj,  can  &e 
expressed  In  the  following  way: 

1  J  ,  4 

m(c,  c)  =  E  E  C,,  •  5  •  cJ  (15) 

1=0  j=0 

Equation  (15)  can  then  be  expressed  In  any  one  of  the  following  forms: 

•  As  a  Surface  Source-Sink  and  Doublet 

After  arbitrarily  choosing  I  =  5  and  J  =  3,  the  strength  of  a 
singularity  distribution  of  surface  source-sink  or  doublet  can  be  expressed 
as 

5  3  ,  . 

m(c,  c)  =  E  E  C..  •  s'  •  (16) 

s  1=0  j=0 

For  the  purpose  of  controlling  the  hull  geometry.  It  is  convenient  to 
rewrite  Equation  (16)  as  follows: 

3  5  ,  i 

m.U,  c)  -  I  (  i  •  e1)  •  (17) 

s  j=0  i=0  1J 

3  j 
■  z  e,  •  r 

j=0  J 

5  . 

with  e.  =  z  (\.  •  £*  (18) 

J  1=0  1J 

where  j  =  0,  1,  2  or  3 

The  singularity  with  eg  of  Equation  (18)  Is  Independent  of  the  depth. 
It  contributes  the  most  to  ship  displacement.  The  singularities  with  e-| 
through  e^  are  more  Influential  on  the  lower  waterlines  and  thus  primarily 
affect  the  shape  of  the  sections. 
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•  As  a  Vertical  Line  Source-Sink  and  Doublet 

The  strength  of  a  vertical  line  source-sink  or  doublet  can  be 
expressed  as 

3  i 

m  (c)  =  z  C.  •  (19) 

v  j=0  J 

t  As  a  Horizontal  Line  Source-Sink  and  Doublet 

The  strength  of  a  horizontal  line  source-sink  or  doublet  can  be 
expressed  as 

5  i 

\U)  =  z  C.  •  c  (20) 

n  i=0  1 

Both  vertical  and  horizontal  line  source-sinks  and  doublets  can  be  placed 
at  any  location  along  the  ship  length  and  depth  that  is  desirable. 

•  As  a  Bottom  Source-Sink  and  Doublet 

(Not  yet  coded  Into  the  computer  program) 

To  obtain  a  flat  bottom,  a  bottom  singularity  distribution  is 
introduced  on  a  horizontal  plane,  nb»  slightly  above  the  required  flat- 
bottom  location.  The  density  distribution  Is  determined  so  that  the  down¬ 
ward  induced  velocity  on  the  flat  bottom  Is  zero.  The  strength  of  a  bottom 
singularity  distribution  may  be  expressed  in  the  following  way: 

5  3  i  , 

mU,  n)  =  z  z  Cs,  •  K  •  nJ  (21) 

1=0  j=0 

To  summarize,  not  only  the  point  source-sink  and  doublet  distributions 
but  also  the  surface  and  line  source-sink  and  doublet  distributions  are 
available  for  use  In  any  combination  with  the  particular  n-surface  chosen. 
Thus,  a  large  family  of  practical  hull  forms  can  be  represented. 

PARAMETRIC  REPRESENTATION  OF  SINGULARITY  DISTRIBUTIONS 

To  improve  control  over  the  hull  geometry,  six  parameters— E,  B,  V, 

Y,  TE,  and  TM— are  chosen  for  each  half  body  (Figure  3). 


Figure  3  -  Singularity  Distribution  on  Forebody 
and  Aftbody 


In  Figure  3,  is  defined  as  the  area  under  the  curve  ff(s)  between 
the  abscissas  a^  and  bf,  (and,  similarly,  Bg  for  the  curve  fg (5 ) ) -  and 
Vg  are  the  moments  of  the  areas  B^  and  Bg,  respectively,  about  the  £  ■  0 
axis.  Ef  and  are  function  values  of  f^(c)  at  £  =  bf,  and  5  =  af, 
respectively.  TEf  and  TMf  are  the  slopes  at  points  b^  and  a^. 

For  surface  source-sink  and  horizontal -line  source-sink,  E  and  TE  are 
grossly  related  to  the  entrance  angle  of  the  waterline;  B  to  the  midship 
section  area;  V  to  the  displacement  volume;  and  Y  and  TM  to  the  parallel 
middle  body.  For  the  surface  doublet  and  horizontal -line  doublet,  V  is 
related  to  the  center  of  buoyancy,  and  B  to  the  displacement  volume. 
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These  six  parameters  are  related  to  the  coefficients  of  the  forebody 
(Equation  (15))  In  the  following  way: 


ej(bf> 


1 


E  Cj  j  •  b» 
1*0  13  f 


<Ef>j 


je}M  •  dt 


fbf 

I  5  •  e^e)  •  ds 


5 

z 

1=0 

5 


z  j  |  i  Cjj  •  (b^  ■  ^ )  =  (B*) 


f'j 


t4t cu  •  <4+1  -  4+2>  -  <v*> 


f'j 


ej(af)  * 

5  i 

E  Ci j  •  a- 

1=0  13  T 

■  <Vj 

se^U) 

5  1-1 

e  1  •  C..  .  bl  1 

1=0  13  T 

■  <TEf>j 

n 

5=bf 

aejU) 

35 

5=af 

5  i  _i 

E  1  •  C..  •  a I 

1=0  13  T 

1  (™f)j 

where  j  *  0,  1,  2,  or  3 


(22 


Equation  (22)  can  then  be  written  In  a  matrix  form  as  shown  In 
Equation  (23).  For  ease  of  referencing,  the  left-hand  portion  of  the 
matrix  representation  Is  noted  as  Matrix  (23a);  the  middle  portion  as 
Matrix  (23b);  and  the  right-hand  portion  as  Matrix  (23c). 


(23a)  (23b)  (23c) 

[P]j  =  [Q]  •  [C]j  (23] 


where  j 


0,  1,2,  or  3, 


The  elements  q^  of  [Q]  are: 

*01  =  bf 

1  fhi+1  J+h 

qli  ‘  TTT  (bf  ‘  af  ' 

„  -  1  /J+2  _  1+2  V 

q21  ‘  r+T  (bf  “  af  ' 

<>31  *  af 
"41  ■  1  •  bf"1 

"51  =  *  '  4'’ 

with  1  -  0,  1,  2,  3,  4,  and  5 

Since  the  determinant  of  [Q]  will  not  be  zero,  [Q]  can  be  inverted 
to  [Q]"1.  For  later  referencing,  [Q]’^  will  be  referred  to  as  Matrix  [24]. 
We  can  write 

[c]j  -  cor1  •  [p]j  (25) 

The  expressions  of  Equations  (22)  through  (25)  apply  to  the  forebody. 
The  aftbody  can  be  expressed  in  a  similar  manner. 
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WAVE  AMPLITUDE  FUNCTIONS 

With  density  distributions  known  and  the  n-surface  given,  the  wave 
amplitude  functions  Ac(e,  Z*  c)  and  As(e,  z%  z)  can  be  obtained  by  a  sur¬ 
face  Integration  over  the  singularity  distribution  area.  The  n-surface 
Is  symmetrical  with  respect  to  n  =  0.  One  half  of  the  distribution  Is  on 
the  positive  side  of  the  n-surface  and  the  other  half  Is  on  the  negative 
side.  For  the  positive  side,  using  Equation  (7),  we  have 


Ac(e,  5.  e) 


1 


M  •  K 


7  *  UTV 


0  •  sec^e  •  e  0 


k  •  ?  •  sec  0 


cos  (Kq  •  z  •  sec  0 


+  KQ  •  n  •  tane  •  sece) 


(26) 


1  7 

=  w  •  m  •  —  •  sec  e  •  e 

2  n 


-k  •  s  •  sec  e 

•  cos  (F^  •  Z  +  F2  •  n) 


where  F,  =  K  •  sec0 
I  0 


F«  *  K  •  sec0  •  tano 
i  0 

For  the  negative  side,  using  Equation  (7),  we  have 


As(0,  z»  c)  *  \  •  m 


K  -  -k  •  c  •  sec  0 
.  sec3e  •  e  0 

IT 


cos  (F,  •  Z  -  F 


2  *  n 
(27) 


) 


Combining  the  equations  of  (26)  and  (27),  we  get  the  amplitude  func¬ 
tion  of  a  cosine  elementary  wave  system  ranging  from  — it/2  to  ir/2. 


i  Kn  -a  -k.  •  z  •  secS 

Ac(0,  z,  z)  =  2  ’  m  *  *  *  secJe  ’  e 


cos  (Fi  •  Z  +  F9  •  n) 


1  ko  3  "ko  ’  5  ’  sec  0 

+  7  ’  m  *  ~n  *  sec  9  *  e  •  cos  (F 


1 


2 

5  -  F 


1 


n) 


i  •  m  •  —  •  sec3©  •  e  0 


-k  •  c  •  sec  0 


[cos  (F]  •  Z  +  F2 


•  r 


+  cos  (F,  •  Z  -  F«  •  n)3 


o  _ 3„  _  "o 


-k  •  c  •  sec  0 


=  m  •  —  •  sec  0  •  e 
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cos  (F1  •  c)  •  cos  (F2  •  n) 
(28) 


t*. 


After  Integration  over  the  distribution  surface,  we  have 


Ac(e) 


I  J 
z  z  C 
i=0  j=0 

rb 

•'A 


..A.  sec3e  r 

^  11  J-lr\ 


0  -k  •  c 

j  -  ' 0  c 


sec2e 


C  .  e 


cos 


(F-j  •  0  •  cos  (F2  •  n)  •  d£ 


K  -  I  J  4  4 

=.  -0  .  sec30  -I  £  C„  •  ZJ  •  xl 
11  i=0  j-0  1J 


dc 


(29) 


where 


Tn 

b 


■k  •  c  •  sec  e 
0  •  dc 


■I. 

=  J  e1  •  cos  (F1  •  0  •  COS  (F2  •  n)  • 


(30) 


(31) 


Likewise,  for  the  sine  amplitude  functions  over  the  range  of  -ir/2  to  n/2, 
we  have 

2 

-k  •  t  •  sec  0 

A  (e)  =  m  •  -  •  sec30  •  e  0  -sin  (F]  •  0  •  cos  (Fg  •  n) 

5 '  '  7T 


k  q  I  J  44 

■A.  .  sec30  •  z  i  C.4  •  ZJ  •  X 
*  i=0  j=0  1J 


where  X, 


■f 


S1  •  sin  (Fi  •  c)  •  cos  (F?  •  n) 


df- 


(32) 


(33) 
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The  representations  of  the  wave  amplitudes  for  the  various  singularity 
distributions  are  derived  In  the  following  ways: 

•  Surface  Source-Sink 

By  substituting  Equation  (13)  into  Equation  (31),  we  have  for  the 
kth  line  segment  (see  Figure  2) 

ocbss  =  f  k+1  s1  •  cos  (F1  .  O-cos  {f2  •  [Bn  •  (A,  •  £  +  A2)  +  d]|< 
fCk+l  , 

=  1  O  •  cos  (F,  •  £)  •  cos  (H  •  £  +  6)  •  d£  (34) 


where  H  *  F«  •  B  •  A, 

2  n  1 

G  =  F2  *  (Bn  *  A2  + 

d  *  distance  from  centerplane  to  c  -  £  plane 


-  }  |  k+1  51  •  jcos  [(H  +  F,) .  {  +  B]  +  cos  [(H  -  F, )  •  c  +  G]  dc 

i  i  r?k+i  < 

=  j  jcos  G  £1  •  cos  [(H  +  F-j)  •  £]  •  d£ 


-  sin  G  • 


+  cos 


-  sin  G  • 


where  1  =  0,  1,  2,  3,  4,  or  5 


k+1  1 

G  •  I  £  •  sin  [(H  +  F, )  •  £]  •  d£ 

J*  k 

ft-k+l  . 

G  •  I  £^  •  cos  C(H  -  F, )  •  £]  •  d£ 

Jh 

fw  i  ) 

G  J  £ '  •  sin  [(H  -  F] )  •  £]  •  d£j 
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Likewise,  Equation  (33)  can  be  expanded  Into 


<xl>» 


c  r5k+i  , 

jcos  6  »j  £  •  sin  [(H  +  F-j)  •  5]  •  d? 
fEk+l  , 

+  sin  G  •  I  5*  cos  C(H  +  F, )  •  c]  •  d£ 

■>h 

fEk+l  , 

cos  G  •  51  •  sin  [(H  -  F, )  •  5]  •  d£ 

J«k 

fEk+l  .  , 

sin  6  “  J  51  •  cos  [(H  -  F^)  •  c]  •  de| 


cos  G  • 


where  1  =  0,  1 ,  2,  3,  4,  or  5 

Solutions  of  Equations  (35)  and  (36)  can  be  obtained  functionally.  By 
summing  the  individual  contributions  of  each  line  segment,  the  total  value 
Is  obtained. 

The  amplitude  functions  of  the  cosine  wave  and  sine  wave  follow: 


5  3 


<Ac<6»ss  '  r  ■  sec°e  ■  JE0  c1j  •  zJ  •  <xc>ss  <3; 

<As<°»ss  ■  i  •  sec3e  •  ,ro  lg  Cfj  •  Z3  •  (*J)M  (3! 

•  Surface  Doublet 

By  partially  differentiating  Equation  (29)  with  respect  to  %  and  n, 
the  amplitude  function  of  a  cosine  wave  can  be  obtained. 

a(A (e))  a(A  (e)) 

<Ac<9»sd  ■  — If-  • xo  +  — 5 —  • vo 


ETO'WT  '-mm!' .  T1JW  -^W  ^  T5  — .  "W 


'."r'tvc7rr' '•«■•■■«-.  x-r'-yr.  • 


k  ,  I  J  ,  3x!  axl 

=  —  •  sec’5©  •  E  E  C< .  •  ZJ  •  (-r=£  •  XD  +  r-2-  .  YD) 
*  1=0  J-0  ij  95  3n 


k  ,  I  J 
—  •  sec  0  •  E  E  C.. 


1=0  j=0 


1j 


(x!) 


c'sd 


(39) 


j  3X  on 

where  <x<^sd  =  af  *  X0  +  3T 


1 


3X 


1 


YD 


•  ”/.' 


[-F,  •  C1 


sin  (F^  •  c)  +  1  •  •  cos  (F^  •  5)] 


cos  (Fg  •  n)  •  d£ 
.b 


YD  J  F2  •  51  •  cos  (F1  •  0  •  sin  (F?  •  n)  •  d£  (40) 


where  1  =  0,  1,  2,  3,  4,  or  5. 

XD  and  YD  are  the  doublet  components  In  the  5-  and  n-dlrectlons, 
respectively,  and  s1  and  are  the  doublet  strengths.  When  the  direction 
of  the  doublet  is  pointing  to  the  ^-direction,  we  have 


XD  =  1.0 
YD  =  0.0 


Likewise,  the  amplitude  function  of  the  sine  wave  is 
9(A  (0))  3(A  (0)) 

<Ve»sd  =  -lr--XD  +  — 1^* YD 

k  n  I  J  ,4 

=  .  sec30  •  E  E  C..  •  ZJ  •  (XM  ,  (41) 

If  j _g  j _q  'J  S  5Q 
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where  (X^)sd  =  XDj  [Fi  •  .  cos  (Fi  •  c)  +  1  •  £  •  sin  (F]  • 


COS  (F«  •  r\)  •  d$ 


Fg  •  •  sin  (F^*  £)  •  sin  (F2  •  n)  •  d£  (42) 


•  Point  Source-Sink 

The  wave  amplitude  functions  of  cosine  and  sine  waves  created  by 
a  traveling  point  source-sink  can  be  deduced  from  Equations  (28)  and  (32), 
respectively. 


k  ,  -k  •  5 • sec  6 

(A  (e))  =  —  •  sec3e  •  m  •  e  0 

c  ps  IT 


cos  (F1  •  5)  •  cos  (F2  •  n)  (43) 


(A ( 0 ) )  n_  -  J 

s  ps  IT 


-  ~k  *  £  *  sec  0 

-  •  sec’5©  •  m  •  e  •  sin  (F^  •  £)  •  cos  (Fg  •  n)  (44) 


•  Point  Doublet 

By  partially  differentiating  Equations  (43)  and  (44)  with  respect 
to  £,  the  wave  amplitude  function  will  be  given  as 

2 

k  >  -k  •  c  •  sec  0 

( Ac ( 6 ) ) pd  =  *  sec  6  •  m  •  e  0  [-F]  .  sin  (F]  .  s) 


•  cos  (F2  •  n)3 


(A_(e))  ,  =  —  •  sec3e  •  m  •  e 


-k  •  c  •  sec  0 


fF,  •  cos  (F,  •  e) 


•  Vertical  Line  Source-Sink 

From  Equations  (29)  and  (32),  the  amplitude  functions  of  cosine 
and  sine  waves  created  by  a  vertical  line  source-sink  may  be  written 


(Me)Lc  *  -r  *  sec3e  • 


c'  "vs 


i= o  J  J> 


-k  •  ?  •  sec  0 
e  0  •  dc 


cos  (F,  •  z)  •  cos  (F,  •  n) 


(47) 


(As(9»vs 


^  •  sec3e  • 

IT 


;  vr\« 

j=0  J  / 


-kQ  •  c  •  sec  0 


dc 


•  sin  (F-j  •  c)  •  cos  (F2  •  n) 


(48) 


•  Vertical  Line  Doublet 

By  partially  differentiating  Equations  (47)  and  (48)  with  respect 
to  £»  the  amplitude  functions  of  cosine  and  sine  wave  created  by  a  vertical 
line  doublet  can  be  expressed 


<Ac<9»vd 


k 


•  sec  0  • 


J=0  J  J, 


•k  •  5  •  sec  0 


dc 


•  [-F1  •  sin  (F1  •  z)  •  cos  (F?  •  n)] 

b 


(49) 
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IT 
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(50) 
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•  Horizontal  Line  Source-Sink 

From  Equations  (29)  and  (32),  the  amplitude  functions  of  cosine 
and  sine  waves  created  by  a  horizontal  line  source-sink  may  be  written 


<\<9»hs 
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(51) 


k  -  5  -k  •  ?  •  sec  e  f  . 

<V9»hs  =  V  •  s“  9  •  C1  •  e  •  5 

1  u  •'a 


•  sin  (F1  •  z)  •  cos  (F2  •  n)  •  d£ 


(52) 


•  Horizontal  Line  Doublet 

By  partially  differentiating  Equations  (51)  and  (52)  with  respect 
to  z,  the  amplitude  functions  of  cosine  and  sine  wave  created  by  a  hori¬ 
zontal  line  doublet  can  be  written 
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COEFFICIENT  OF  WAVEMAKING  RESISTANCE 

Equation  (8)  Indicates  the  wavemaking  resistance  of  a  ship.  A  non- 
dimensional  coefficient,  C  ,  of  wavemaking  resistance  Is  defined  as 

rt 

2  - 

c“  '  777  T?  =  s  f 2  (Ac(s)2  +  as(9)2>  • sec3e  ' de  (55) 

1/2  •  p  •  U  •  L  J q 

where  L  Is  a  suitably  chosen  length. 

Since  the  free-surface  disturbance  Is  assumed  to  be  small  where  the 
boundary  condition  Is  linearized  In  the  theory,  the  total  wave  amplitudes 
A  (e)  and  A  (e)  are  the  sums  of  those  various  singularity  elements 

C  j 

Ac(e)  =  (Ac(e))ss  ♦  (Ac(e))sd  +  <Ac(e))ps  ♦  (Ac(e))pd  ♦  (Ac(o))vs 

+  (Ac(e))vd  ♦  (Ac(e))hs  ♦  (Ac(e))M  (56) 

As(e)  -  (As(e))ss  ♦  <As(e))sd  +  (As(e))ps  +  (As(e))pd  +  (As(e))vs 

+  (As(e))vd  *  (As(e))hs  +  (As(e))hd  (57) 

OPTIMIZATION  OF  WAVEMAKING  RESISTANCE 

Calculating  wavemaking  resistance  from  a  set  of  given  singularity 
distributions  has  already  been  discussed  In  great  detail.  The  Plen  pro¬ 
gram's  main  objective  is  to  obtain  a  set  of  singularity  distributions 
that  will  generate  a  hull  geometry  which  satisfies  the  design  conditions 
and  at  the  same  time  has  good  resistance  performance. 

Two  approaches  to  optimization  are  possible.  In  one  of  these,  a  hull 
form  Is  already  known,  and  the  singularity  distribution  representing  the 
hull  form  can  be  modified.  In  this  case,  the  objective  is  to  improve  the 
resistance  performance  in  the  following  way,  with  or  without  specifying 
changes  In  the  hull  beam  and  displacement.  The  modifying  singularity  dis¬ 
tribution  (Equation  (15))  is  introduced,  using  only  one  value  of  j  at  a 
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time.  The  coefficient  In  the  chosen  modifying  singularity  distribution 
Is  determined  by  optimizing  the  wavemaking  resistance  of  all  of  the  sin¬ 
gularity  distributions  within  the  constraints  Imposed  on  the  modifying 
singularity  distribution.  When  this  has  been  done,  the  hull  geometry  is 
computed  to  determine  whether  or  not  It  is  practical.  If  it  is,  another 
modifying  singularity  distribution  with  a  different  j  value  is  chosen  and 
the  optimization  computation  repeated.  If  It  is  not,  a  different  set  of 
constraints,  or  even  a  different  modifying  singularity  distribution,  is 
tried.  If  neither  of  these  techniques  provides  a  more  satisfactory  result, 
the  result  obtained  before  the  current  optimization  attempt  should  be  used. 

The  procedure  just  described  derives  a  new  design  by  improving  an 
existing  design  through  optimization.  This  method  works  well  when  the 
parent  form  to  be  used  has  a  known  singularity  distribution.  The  new 
design  is  achieved  by  Imposing  appropriate  constraints  on  the  modifying 
singularity  distributions  of  the  parent  form.  When  there  is  no  appropriate 
parent  form  available,  the  entire  singularity  distribution  for  the  new 
design  has  to  be  obtained  from  the  optimization  process.  It  is  advisable 
to  start  the  first  step  of  the  optimization  procedure  with  a  smaller  dis¬ 
placement  and  beam  than  would  likely  be  suitable.  The  modifying  singu¬ 
larity  distribution  chosen  for  subsequent  steps  of  the  optimization  and 
the  constraints  specified  will  then  be  guided  by  the  result  of  the  pre¬ 
vious  step.  The  optimization  process  used  is  as  follows: 


For  a  value  of  j  corresponding  to  each  term  of  Equation  (15),  there  is  a 
wave  amplitude  function  d.  •  •  ZJ  •  X^.  By  applying  Equations  (30)  and 

*  J 

(31)  an  arbitrary  wave  amplitude  function  defined  as 


is  introduced.  This  wave  system  is  designed  to  partially  eliminate  the 
wave  system  generated  by  the  existing  wave  amplitudes  of  the  original 
singularity  distributions  (Equations  (56)  and  (57)). 
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The  six  coefficients  of  Equation  (58) 


’ij 


1  =  0,  1,  2,  3,  4,  5 


(59) 


where  j  =  0,  1,  2,  or  3,  are  to  be  selected  under  a  set  of  design  con¬ 
straints  and  by  an  optimization  procedure. 

Assume  from  Equations  (55)  and  (58) 

o  £ 
i.c  r  7T 
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c  p  2  g  g 

M  [(  E  d.^  •  Zj  .  +  A  (e))2  +  (  Z  d.,  •  Z 1  •  xj  +  A  (0))2] 

tt  J  i=Q  c  c  i=0  1J  s  s 


•  sec  e  •  de 
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+  Ac(e)[  de 
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+  27 


if 


5  5 

t  [  £  (di.j  •  d..) 

k=0  1=0  KJ 


*1  •  ‘s +  2  •  dkj  •  ‘s  • 


+  As(0)[  d0 


(60) 


where  Ag  and  A£  are  free-wave  amplitudes  of  the  original  singularity  dis¬ 
tribution  and 


I. 


•  sec**  •  de 


TT 

2 


=  I  Z^  •  Xg  •  sec"*  •  de, 
J  n 


(61) 


where  j  =  0,  1,  2,  or  3 
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The  variables  Z^,  and  are  defined  in  Equations  (30),  (31),  and 

V  d 

(33),  respectively. 

After  the  terms  shown  in  the  first  integral  of  the  expanded  form  of 
Equation  (60)  are  expanded,  they  can  be  rearranged  as  a  7X7  matrix 


where  j  =  0,  1,  2,  or  3  and  the  elements  b.  .  of  [B].  are 


b01  -  d0j  •  ‘c  <d1j  •  +  Ac) 

bH  *  dlj  •  <d1j  •  ‘c  +  Ac> 

b21  ■  d2j  •  ‘c  <d1j  '  ‘c  +  Ac) 

b3i  ■  d3j  •  i  <dij  •  *;  -  ac> 

b4i  *  d4j  •  ‘c  (d1j  •  ‘c  +  Ac> 

be*  *  dc<  *  .  *  ^  *  +  A  ) 


Ac  <di: 


‘c  +  Ac> 


where  i  =  0,  1,  2,  3,  4,  and  5. 

After  separating  the  dfc.  from  each  term,  the  remaining  terms  can  be 
written  as 


t°  (t®  +  Ac)  t®  (tj  *  Ac) 
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In  a  similar  way,  the  terms  of  the  second  integral  of  Equation  (60) 
are  expanded  and  then  rearranged  as  a  7X7  matrix  form 
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The  elements  b^  of  [B]j  are 
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where  i  =  0,  1,  2,  3,  4,  and  5. 

After  separating  the  d^  from  each  term,  the  remaining  terms  can  be 
written  as 


‘s  <‘s  +  As> 
‘1  »*J  +  As> 
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(65) 


The  six  coefficients  defined  In  Equation  (58)  are  the  only  unknowns 
of  the  Equation  (60).  These  six  numbers  can  be  expressed  parametrically 
as  E,  B,  V,  Y,  TE,  TM  by  Equation  (22).  An  exponential  random  search 
method^  was  used  to  search  for  a  minimum  value  of  by  randomly  varying 
the  six  parameters  within  the  range  of  their  search  limits  (or  design  con¬ 
straints). 


34 


Equation  (58),  the  added  wave  amplltuue  function,  may  be  any  one  of 
the  following  singularity  distributions  In  which  j  will  have  the  value 
0,  1,  2,  or  3: 

•  surface  source-sink 

•  surface  doublet 

•  horizontal  line  source-sink 

•  horizontal  line  doublet 

Singularity  distribution  of  vertical  line  and  bottom  can  easily  be 
Included  In  the  optimization  process  at  some  future  time. 

After  one  optimization  Iteration,  values  of  the  six  numbers  (d^) 

5 

just  selected  are  lumped  Into  E  C..  of  Ar  and  Ac  of  Equation  (60)  to 

1=0  1J  u  s 

form  the  new  singularity  distribution.  The  second  optimization  iteration 
Is  then  begun,  with  j  =  0,  1,  2,  or  3,  to  find  another  set  of  d...  This 

*  J 

process  will  be  repeated  over  and  over  until  the  value  of  Cw  can  no 
longer  be  reduced.  Figure  4  shows  the  flow  diagram  of  this  optimization 
process. 
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*Singularity  distributions  available: 

Surface  source-sink,  forebody 
Surface  Source-sink,  aftbody 
Surface  doublet,  forebody 
Surface  doublet,  aftbody 
Horizontal  line  source-sink,  forebody 
Horizontal  line  source-sink,  aftbody 
Horizontal  line  doublet,  forebody 
Horizontal  line  doublet,  aftbody 


Figure  4  -  Flow  Chart  of  the  Optimization  Process 
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STREAMLINE  FUNCTIONS  AND  VELOCITY  COMPONENTS 


The  velocity  potential  $  is  given  by  the  three  dimensional  singular¬ 
ity  distribution  regarding  the  distribution  as  symmetrical  to  the  free 
surface.  The  velocity  components  u,  v,  and  w  are  induced  by  the  singu¬ 
larity  distributions  at  any  field  point  (x,  y,  z). 


•  Surface  Source  Sink 
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where  r  =  ^  (x-c)2  +  (y-n)2  +  (z-c)‘ 
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v  _  .  Z£d  .  J_  .  z*LL 
pd  3  y  4ir 


w  -  8*Pd  -  JL  •  f-3m  •  (x-e)  ‘  (z-c)l 
"pd  '  '  11  "  4.  [  r5  J 

where  m  is  the  strength  of  the  doublet  divided  by  ship  speed. 


•  Vertical  Line  Source  Sink 


where  the  line  source- sink  extends  in  depth  from  a  to  b  and  is  located 
perpendicular  to  the  E-n  plane. 
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where  the  line  doublet  extends  In  depth  from  a  to  b  and  is  located 
perpendicular  to  the  z-n  plane. 


•  Horizontal  Line  Source-Sink 
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where  the  horizontal  line  source-sink  extends  in  the  x-direction  from  a 
to  b  and  is  located  perpendicular  to  the  n-5  plane. 


e  Horizontal  Line  Doublet 
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where  the  line  doublet  extends  in  the  x-di recti  on  from  a  to  b  and  is 
located  perpendicular  to  the  n-€  plane. 

STREAMLINE  TRACING 

By  introducing  the  velocity  components  derived  in  the  previous 
section  into  Equation  (10)  and  then  integrating,  we  can  obtain  the 
streamlines  of  the  flow.  Equation  (10)  can  be  solved  numerically  by 
Range-Kutta  methods  when  the  starting  points  of  the  streamlines  are 
known.  To  obtain  a  desired  accuracy  of  the  streamlines,  however,  the 
Kutta-Merson  method  of  integrating  the  differential  equations  has  been 
chosen.  This  method  automatically  regulates  the  step-length  used  in 
tracing  streamlines  so  that  small  step-lengths  are  used  in  regions  where 
there  are  extreme  changes  in  curvature— such  as  close  to  the  bow  contour, 
for  instance— and  larger  step-lengths  are  used  in  regions  where  the 
changes  in  u,  v,  and  w  are  gradual. 

To  determine  proper  starting  points  for  the  streamlines,  it  has 
been  foum.1  advantageous  to  determine  the  stagnation  point  on  the  bow 
contour  and  then  to  trace  the  bow  contour  streamline  in  the  center  plane 
above  and  below  the  stagnation  point.  Starting  points  for  the  stream¬ 
lines  can  afterwards  be  chosen  on  the  bow  contour,  slightly  offset  from 
the  centerplane.  The  forebody  is  obtained  by  choosing  an  appropriate 
number  of  starting  points  and  then  tracing  streamlines  from  the  bow 
contour  rearward. 


IV.  DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


PROGRAM  STRUCTURE 


Figure  5  shows  the  program  structured  as  overlays.  The  program 
reads  Input  data  with  overlay  (1,0)  In  core.  Then  overlay  (2,0)  Is 
brought  Into  core  to  compute  the  wavemaking  resistance  from  Input 
singularity  distributions.  The  wavemaking  resistance  may  be  optimized 
to  obtain  the  new  singularity  distributions  under  a  set  of  design  con¬ 
straints.  At  this  stage,  overlay  (3,0)  may  be  brought  into  core  to 
replace  overlay  (2,0).  Streamlines  corresponding  to  the  input  or 
optimized  singularity  distributions  can  now  be  traced  to  see  whether  the 
ship  hull  form  they  represent  is  practical. 

(0,0) 


(1,0) 


PIEN  -  Main 
Control  Program 


_ I _ •  '  ' 

PCIPT  -  Reads 

PCWMR  -  Computes 

_ i _ " 

PCSTL  -  Traces 

and  Prints 

and  Optimizes  Wave- 

Ship  Hull  Form 

Input  Data 

making  Resistance 

i 


HWWMR  -  Computes 
Wavemaking  Resistance! 


HWOPT  -  Optimizes 
Wavemaking  Resistance 


Figure  5  -  Program  Overlay  Structure 


Figures  6  through  9  which  follow  show  the  organization  of  programs 
within  each  overlay. 


I 
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* 


*  NSRDC  Library  Subroutines. 

Figure  6  -  Block  Diagram  For  (1,0)  Overlay  - 
Input  Data  Processing 
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♦Subroutines  described  by  Gray.^ 


Figure  8  -  Block  Diagram  for  (2,0)  Overlay  - 
Wavemaking  Resistance  Optimization 


SUBROUTINE  DESCRIPTIONS 


The  following  pages  contain  detailed  functional  descriptions  of  the 
various  subroutines  of  the  Plen  Wavemaking  Resistance  Computation  Program. 
A  complete  alphabetical  listing  of  the  subroutines  by  overlay  follows: 
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*  Wavemaking  Resistance  Computation 
♦♦Wavemaking  Resistance  Optimization 

Preceding  page  blank 


PROGRAM  PIEN 


Function: 

This  program  Is  the  main  program  In  the  (0,0)  overlay. 


Calling  Format: 

PROGRAM  PIEN  (OUTPUT,  TAPE6=0UTPUT,  INPUT,  TAPE5=INPUT) 


Detailed  Description: 

Program  PIEN  first  Initializes  several  arrays  and  then  calls  In  the 
programs  in  the  (1,0)  overlay. 


Subroutines  Called: 
PCIPT 


SUBROUTINE  HWCTP 


Function: 

This  subroutine  computes  solutions  of  the  six  linear  simultaneous 
equations. 

Calling  Format:  +  *  i  I  +  + 

CALL  HWCTP  (L,  SV,  KP,  D,  R,  PV) 

Description  of  Parameters: 

L  An  Indicator: 

=  1  The  forebody 
=  2  The  aftbody 

SV(6)  Contains  Matrix  (23C) 

KP  A  printout  indicator: 

t  1  No  printout 
=  1  Print  array  PV 

D,R  D=a,  R=b  (Figure  3) 

PV (6)  Contains  Matrix  (23a) 

Detailed  Description: 

Subroutine  HWCTP  solves  Equation  (23).  If  the  value  in  KP  is  1, 
the  program  prints  out  the  values  of  Matrix  (23a). 

Subroutines  Called: 

None 
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SUBROUTINE  HWHPI 


Function: 

This  routine  reads  input  data  in  "point"  form  for  singularity 
distribution  of  surface  source-sink  and  doublets. 


Calling  Format:  +  f 

CALL  HWHPI  (CLH,  SCH,  L) 


Description  of  Parameters: 


CLH (10, 4)  Refer  to  the  description  of  CLSH  or  CLDH  in  Appendix  A. 
SCH(2)  Scaling  factors 

L  An  indicator: 

=  1  The  forebody 
=  2  The  aftbody 


Detailed  Description: 

Subroutine  HWHPI  reads  in  the  strength  of  a  singularity  distribution 
in  "point"  form  and  then  calls  the  POINTS  subroutine  to  compute  the 
coefficients  of  a  5th-degree  polynomial  which  is  an  approximation  based 
upon  the  "point"  data. 

Subroutine  Called: 

POINTS 
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SUBROUTINE  HWLHI 
Function: 

This  subroutine  reads  and  prints  the  input  data  for  horizontal 
source-sink  and  doublets. 

Calling  Format: 

CALL  HWLHI 

Detailed  Description: 

Subroutine  HWLHI  reads  in  the  Integration  interval  for  the  singular¬ 
ity  distribution,  and  then  calls  the  HWMTX  subroutine  to  compute  Matrix 
(24)  for  both  the  fore  and  aftbody  section.  It  then  reads  the  strength 
of  singularity  distribution  presented  in  one  of  three  forms:  as 
coefficients,  as  parameters,  or  as  "points".  If  in  parameter  form, 
subroutine  HWPLI  is  called;  if  in  "point"  form,  subroutine  HWHPI  is 
called.  The  program  writes  out  all  of  the  input  data  read. 

If  11<KWR<20  or  if  31<KWR<40,  the  program  prints  only  the  input 
data  contained  in  "COMMON"  storage. 

Subroutines  Called: 


HWCTP,  HWMTX,  HWHPI,  HWPLI,  HWPTS 


SUBROUTINE  HWLVI 


Function; 

This  subroutine  reads  and  prints  Input  data  for  vertical  line 
source-sinks  and  doublets. 

Calling  Format; 

CALL  HWLVI 

Detailed  Description; 

Subroutine  HWLVI  reads  In  and  prints  out  the  CLSV  and  CLDV  arrays 
which  contain  Information  pertaining  to  the  vertical  line  source-sinks 
and  doublets. 

Subroutines  Called; 

None 


SUBROUTINE  HWMTX 


Function: 

Sets  up  six  simultaneous  linear  equations  and  computes  the  inverse 
matrix  for  the  coefficients  of  the  equations. 


Calling  Format:  +  +  +  4. 

CALL  HWMTX  (L,  CMS,  D,  R) 

Description  of  Parameters: 

L  The  index  of  the  three-dimensional  array  CMS,  (1<KL<4) 

CMS(6,6,4)  Refer  to  the  description  of  CMH  or  CMS  in  Appendix  A. 
D,R  D=a,  R=b  (See  Figure  3) 

Detailed  Description: 

HWMTX  first  fills  the  CMS  array  with  Matrix  (23b)  for  surface 
source-sink/doublet  or  horizontal  line  source-sink/doublet,  then  calls 
Subroutine  MATINS  to  compute  Matrix  (24). 

Subroutines  Called: 

MATINS 
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SUBROUTINE  HWPLI 


Function: 

This  subroutine  reads  Input  data  In  the  parametric  form  for  the 
horizontal  source-sinks  or  doublets. 


Calling  Format:  +  +  +  f  +  + 

CALL  HWPLI  (CLH,  CMH,  PV,  SCH,  L,  K) 


Description  of  Parameters: 

CLH(10,4)  Refer  to  description  of  CLSH  or  CLDH  In  Appendix  A. 
CMH(6,6,4)  Refer  to  Appendix  A. 

PV(6,2,4)  Contains  Matrix  (23a)  for  either  the  horizontal  source- 
sink  or  the  doublet. 


SCH(2)  Scaling  factors 

L  An  indicator: 

=  1  The  forebody 
=  2  The  aftbody 


K  An  Indicator: 

»  0  Horizontal  source-sink 
=  2  Horizontal  doublets 


Detailed  Description: 

Subroutine  HWPLI  reads  In  six  parameters  of  Matrix  (23a)  for 
horizontal  source-sink  or  doublet  and  then  calls  HWPTC  subroutine  to 
convert  the  six  parameters  Into  coefficient  form  (Matrix  (23c)). 


Subroutine  Called: 
HWPTC 
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SUBROUTINE  HWPOI 
Function: 

This  routine  reads  Input  data  In  "point"  form  for  singularity 
distribution  of  surface  source-sink  and  doublets. 

Calling  Format:  + 

CALL  HWPOI  (CES,  SCS,  L) 

Description  of  Parameters: 

CES  (I,J,L)  Refer  to  description  of  CESS  or  CESD  In  Appendix  A. 

SCS(4,2)  Scaling  factors 

L  An  indicator: 

=1  The  forebody 
-2  The  aftbody 

Detailed  Description: 

Subroutine  HWPOI  reads  In  the  strength  of  a  singularity  distribution 
In  "point"  form  and  then  calls  the  POINTS  subroutine  to  compute  the 
coefficients  (Matrix  (23c))  of  a  5th-degree  polynomial  which  Is  an 
approximation  based  upon  the  "point"  data. 

Subroutine  Called: 


POINTS 
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SUBROUTINE  HWPRI 
Function: 

This  subroutine  reads  input  data  in  the  parametric  form  for  the 
surface  source-sinks  or  doublets. 

Calling  Format:  f  +  +  +  t  *  + 

CALL  HWPRI  (XS,  CMS,  CES,  PV,  SCS,  L,  K) 

Description  of  Parameters: 

XS(2,2)  Refer  to  the  description  of  XSS  or  XSD  in  Appendix  A. 
CMS(6,6,4)  Refer  to  Appendix  A. 

CES(6,4,4)  Refer  to  description  of  CESS  or  CESD  in  Appendix  A. 

PV(6,2,2)  Contains  Matrix  (23a)  for  either  the  surface  source-sink 
or  the  doublet. 

SCS(4,2)  Scaling  factors 

L  An  indicator: 

=  1  The  forebody 
=  2  The  aftbody 

K  An  indicator: 

=  0  Surface  source-sink 
*  2  Surface  doublets 

Detailed  Description: 

Subroutine  HWPRI  reads  in  six  parameters  of  Matrix  (23a)  for  surface 
source-sink  or  doublet  and  then  calls  HWPTC  subroutine  to  convert  the 
six  parameters  into  coefficient  form  (Matrix  (23c)). 

Subroutine  Called: 

HWPTC 
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SUBROUTINE  HWPTC 


Function: 

This  subroutine  converts  solutions  of  a  set  of  six  simultaneous 
linear  equations  Into  coefficients. 


Calling  Format: 


4-  4-  +  4-  I 


CALL  HWPTC  (L,  PV,  KP,  D,  R,  SV,  J,  CMS) 


Description  of  Parameters: 


L  An  Indicator: 

*  1  The  forebody 

*  2  The  aftbody 


PV (6)  Contains  Matrix  (23a) 


KP  A  printout  direction: 

f  1  No  printout 

=  1  Both  parameters  and  coefficients  of  the  six 
equations  are  to  be  printed 

D,R  D=a,  R=b  (See  Figure  3). 

SV (6 )  Contains  Matrix  (23c) 

J  Contains  the  value  of  j+1  of  Equation  (18) 

CMS  (6,6)  Refer  to  description  of  CMS  or  CMH  in  Appendix  A. 


Detailed  Description: 


This  subroutine  solves  Equation  (25)  using  Matrices  (24)  and  (23a). 


If  the  value  in  KP  is  1,  the  program  prints  out  the  values  of  Matrices 
(23a)  and  (23c). 


Subroutines  Called: 
None 
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SUBROUTINE  HWPTI 
Function: 

This  subroutine  reads  input  data  for  point  source-sink  and  doublet. 

Calling  Format: 

CALL  HWPTI 

Detailed  Description: 

Subroutine  HWPTI  reads  the  scaling  factors  (PSL(I),  1=1,5)  for 
each  group  of  point  source-sinks  or  doublets,  and  scales  each  point 
source-sink  or  doublet  as  follows: 

PINS(1 ,*)  =  PINS(1,*)*PSL(1)+PSL(5) 

PINS(2 ,*)  =  PINS(2,*)*PSL(2) 

PINS(3 ,*)  =  PINS(3,*)*PSL(3) 

PINS(4 ,*)  =  PINS(4,*)*PSL(4) 

It  then  can:  Subroutine  SORTR  to  sort  the  array  PINS  or  PIND  into 
descending  order  according  to  the  5-value  of  each  data  point. 

Subroutine  Called: 


SORTR 


SUBROUTINE  HWPTS 
Function: 

This  routine  Interprets  a  5th-degree  polynomial  equation. 

Calling  Format:  +  +  +  + 

CALL  HWPTS  (L,  SV,  J,  KPL) 

Description  of  Parameters: 


L 

An  Indicator: 

*  1  The  forebody 
=  2  The  aftbody 

SV(6) 

Contains  Matrix  (23c) 

J 

Value  of  j+1  of  Equation  (18) 

KPL 

Not  used 

Detailed  Description: 

Subroutine  HWPTS  evaluates  and  prints  out  a  5th-degree  polynomial 
(Equation  (18))  for  the  Interval  x  -  -1.3  to  x  ■  1.3  with  increment  size 
0.1. 


Subroutines  Called: 
None 


( 
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SUBROUTINE  HWSFI 


Function: 

This  subroutine  reads  and  prints  the  Input  data  for  surface  source- 
sink  and  doublets. 

Calling  Format: 

CALL  HWSFI 

Detailed  Description: 

Subroutine  HWSFI  reads  In  the  variable  defining  the  n-surface  and 
the  Integration  for  the  singularity  distribution,  and  then  calls  the 
HWMTX  subroutine  to  compute  Matrix  (24)  for  both  the  fore  and  aftbody 
sections.  It  then  reads  the  strength  of  singularity  distribution  supplied 
In  one  of  three  forms:  as  coefficient,  as  parameters,  or  as  "points". 

If  the  input  data  Is  in  parameter  form,  the  subroutine  HWRPI  is  called; 

If  in  "point"  form,  the  subroutine  HWPOI  is  called.  When  the  value  In 
NIT  Is  greater  than  one,  the  n-surface  will  be  subdivided  into  NIT 
subsegments  by  calls  to  HWXCD.  The  program  writes  out  all  of  the  input 
data  read. 

If  1 1 <KWR<20  or  If  31<KWR<40,  the  program  prints  only  the  input 
data  contained  In  COMMON  storage. 

Subroutines  Called: 

HWCTP ,  HWMTX,  HWPOI,  HWPRI,  HWPTS,  HWXCD 
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SUBROUTINE  HUXCD 


Function: 

This  subroutine  approximates  an  n-surface  by  using  a  number  of 
straight  line  segments. 

Calling  Format:  +  +  *  +  +  * 

CALL  HUXCD  (XI,  X2,  Kl,  K2,  XD,  L) 

Description  of  Parameters: 

XI, X2  The  beginning  and  the  end  of  a  straight  line  segment 

Kl  The  number  of  the  straight  line  segments  (1<K1<4) 

K2  An  indicator: 

=  1  Surface  source-sink 
=  2  Surface  doublet 

XD  Distance  between  XI  and  X2  (X2-X1) 

L  An  indicator: 

=  1  The  forebody 

=  2  The  aftbody 

Detailed  Description: 

Subroutine  HUXCD  computes  Aj  and  A2  of  Equation  (13)  and  stores  the 
results  into  the  AAA  and  BBB  arrays  respectively.  For  description  of 
arrays  AAA  and  BBB,  refer  to  Appendix  A. 

Subroutines  Called: 

None 


■ 
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SUBROUTINE  PCIPT 


Function: 

This  subroutine,  the  main  controlling  routine  In  the  (1,0)  overlay, 
acts  as  the  executive  program  for  Input-data  processing. 

Calling  Format: 

CALL  AETSKC(5HPCIPT) 

Detailed  Description: 

Subroutine  PCIPT  reads  In  the  printout  titles,  the  ten  calculation 
flags,  the  Froude  numbers,  and  other  variables.  Depending  upon  the 
calculation-flag  values,  it  may  then  call  upon  other  subroutines  to 
read  Input  data  for  a  surface  source-doublet,  a  point  source-doublet, 
a  vertical-line  source-doublet,  or  a  horizontal -line  source-doublet. 

It  will  then  transfer  control  to  PCWMR  In  the  (2,0)  overlay  or  to  PCSTL 
In  the  (3,0)  overlay,  depending  upon  the  value  of  the  Input  parameter  KWR. 

Subroutines  Called: 


HWSFI ,  HWPTI,  HWLVI ,  HWLHI,  PCWMR,  PCSTL 
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SUBROUTINE  SORTR 


Function: 


This  subroutine  sorts  an  array  into  descending  order. 


Calling  Format:  +  +  + 

CALL  SORTR(KP,PIN,NP) 


Description  of  Parameters: 

KP(100)  The  array  PIN  in  its  original  order. 

PIN(4,100)  Refer  to  the  description  of  PINS  of  PIND  in  Appendix  A. 
NP  Refer  to  description  of  NPS  of  NPD  in  Appendix  A. 


Detailed  Description: 

Subroutine  SORTR  arranges  the  contents  of  the  PIN  array  in  a  descending 
order  according  to  the  first  element  PIN( 1 ,*)  of  the  array. 


Subroutines  Called: 


None 
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FUNCTION  FRICT 

Function: 

This  function  computes  the  coefficient  of  frictional  resistance. 

Calling  Format: 

FRICT  (VS,  SL,  I FT) 

Description  of  Parameters: 

VS  Ship  speed 
SL  Ship  length 
IFT  Indicator: 

=1  ATTC  formula  used 
=2  ITTC  formula  used 

Detailed  Description: 

Function  FRICT  computes  the  coefficient  of  frictional  resistance  by 
using  the  formulas  in  Section  3,  Chapter  VII  of  Comstock.7 

Subroutine  called: 

None 


i 
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SUBROUTINE  HWCMR 
Function: 

This  subroutine  computes  the  coefficient  of  wavemaking  resistance. 


Calling  Format: 

CALL  HWCMR 

Detailed  Description: 

Subroutine  HWCMR  first  sums  up  the  various  wave  amplitude  functions. 
After  squaring  Ac  (Equation  (56))  and  Ag  (  Equation  (57)),  it  calls  the 
function  SIMPUN  for  the  integration  of  the  sum  of  these  squares  to 
obtain  the  coefficient  of  wavemaking  resistance  as  shown  in  Equation  (55). 

Subroutine  Called: 

SIMPUN 
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SUBROUTINE  HWLXX 
Function: 

This  routine  computes  the  Integration  term  on  the  right  side  of 
Equations  (51)  and  (52)  or  (53)  and  (54). 


Calling  Format:  +  + 

CALL  HWLXX  (XS,  XE, 


IND,  L) 


Description  of  Parameters: 


XS.XE 

IND 


Starting  and  ending  point  of  the  Integration  Interval 
In  the  5-dlrectlon 

An  Indicator: 


=  1  Surface  source-sink 
*  2  Surface  doublet 


L  An  Indicator: 

=  1  The  forebody 
=  2  The  aftbody 


Subroutine  Called: 


PPP 
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SUBROUTINE  HWPT2 


Function: 

This  subroutine  prints  amplitude  functions  of  the  cosine  and  sine  wave. 

Calling  Format: 

CALL  HWPT2  (KFW) 


Description  of  Parameters: 

KFW  Index  to  the  number  of  different  wave  amplitudes  to  be 
printed. 

Detailed  Description: 

HWPT2  first  prints  values  at  different  degrees  of  the  cosine  wave 
amplitude  functions  of  surface  source-sink  (Eq . (37 ) ) ,  surface  doublet 
(Eq.  (39)),  point  source-sink  (Eq.  (43)),  point  doublet  (Eq.  (45)), 
vertical  line  source-sink  (Eq.  (47)),  vertical  line  doublet  (Eq.  (49)), 
horizontal  line  source-sink  (Eq.  (51))  and  horizontal  line  doublet 
(Eq.  (53))  and  then  prints  values  at  different  degrees  of  the  sine  wave 
amplitude  functions  of  surface  source-sink  (Eq.  (38)),  surface  doublet 
(Eq.  (41)),  point  source-sink  (Eq.  (44)),  point  doublet  (Eq.  (46)),  vertical 
line  source-sink  (Eq.  (48)),  vertical  line  doublet  (Eq.  (50)),  horizontal 
line  source-sink  (Eq.  (52))  and  horizontal  line  doublet  (Eq.  (54)). 

Subroutine  Called: 

None 


* 

* 

i 
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SUBROUTINE  HWPT3 
Function: 

This  routine  prints  the  speed-length  ratio,  coefficients  of 
frictional  and  wavemaking  resistance,  and  horsepower  at  different  ship 
speed. 

Calling  Format: 

CALL  HWPT3 

Detailed  Description: 

Subroutine  HWPT3  first  reads  from  input  the  length  of  ship, 
wetted  surface  of  the  ship,  and  other  information.  It  then  enters  into 
a  DO  loop.  In  the  DO  loop,  HWPT3  prints  the  Froude  number,  speed-length 
ratio,  coefficients  of  frictional  and  wavemaking  resistance,  and 
horsepower  at  different  ship  speeds. 

Subroutine  Called: 

FRICT 
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SUBROUTINE  HWWMR 
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Function: 

This  subroutine  Is  the  main  execution  program  for  the  wavemaking 
resistance  computation. 

Calling  Format: 

CALL  HWWMR 

Detailed  Descriptions: 

Subroutine  HWWMR  contains  two  main  Iterative  loops.  The  outer  loop 
Is  for  the  computation  of  the  wavemaking  resistance  at  specific  Froude 
numbers  beginning  with  FDS  and  ending  with  FDE  with  the  Increment  FDI. 

The  Inner  loop  Is  for  the  computation  of  the  amplitude  function  of  the 
waves  at  the  directional  angles  between  0  and  85  degrees;  the  cosine  and 
sine  components  of  the  wave  amplitude  functions  are  then  computed  by  the 
successively  called  subroutines  HWXSS,  HWXSD,  HWXPS,  HWXPD,  HWXLSV,  HWXLDV, 
HWXLSH,  and  HWXLDH.  These  wave  amplitudes  are  stored  In  the  WAMP  array. 
The  last  row  of  the  7X7  triangular  matrix  TSC,  which  corresponds  to  the 
Matrices  of  Equations  (63)  and  (65),  Is  also  filled  here  for  the  optimi¬ 
zation  of  wavemaking  resistance.  When  the  Inner-loop  computation  Is 
complete,  the  wavemaking  resistance  (Equation  (55))  Is  computed  by  a  call 
to  HWCMR.  The  computations  are  saved  In  the  CRW  array. 

Before  the  program  returns,  the  wavemaking  resistance  and  wave 
amplitude  are  printed  out. 

Subroutines  Called: 

HWXSS,  HWXSD,  HWXPS,  HWXPD,  HWLSV,  HWLDV,  HWLSH,  HWLDH ,  HWCMR, 

HWPT2 ,  HWPT3,  HWZSD 
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SUBROUTINE  HUXLDH 
Function; 

This  routine  computes  the  cosine  and  sine  components  of  wave 
amplitudes  generated  by  the  horizontal  line  doublet. 

Calling  Format; 

CALL  HUXLDH 

Detailed  Description; 

This  routine  computes  Equations  (53)  and  (54)  and  stores  the  computed 
wave  amplitudes  in  the  WAMP  array.  The  first  six  rows  of  the  TSC  array 
for  the  optimization  of  wavemaking  resistance  are  also  filled  if  the 
value  In  KLP  Is  1. 

Subroutines  Called; 

HWLXX 
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SUBROUTINE  HUXLDV 


Function: 

This  routine  computes  the  cosine  and  sine  components  of  wave 
amplitude  generated  by  the  vertical  line  doublet. 

Calling  Format: 

CALL  HWXLDV 

Detailed  Description: 

Subroutine  HWXLDV  computes  Equations  (49)  and  (50)  and  stores  the 
computed  wave  amplitudes  Into  the  WAMP  array. 


Subroutines  Called: 


SUBROUTINE  HWXLSH 
Function; 

This  routine  computes  the  cosine  and  sine  components  of  wave 
amplitudes  generated  by  the  horizontal  line  source-sink. 

Calling  Format: 

CALL  HWXLSH 

Detailed  Description: 

Subroutine  HWXLSH  computes  Equations  (51)  and  (52)  and  stores  the 
computed  wave  amplitudes  In  the  WAMP  array.  The  first  six  rows  of  the 
TSC  array  for  the  optimization  of  wavemaking  resistance  are  also  filled 
If  the  value  In  KLP  Is  1. 

Subroutines  Called: 


HWLXX 


SUBROUTINE  HWXLSV 
Function: 

This  routine  computes  the  cosine  and  sine  components  of  wave 
amplitudes  generated  by  the  vertical  line  source-sink. 


Calling  Format: 


CALL  HWXLSV 


Detailed  Description: 


This  routine  computes  Equations  {47)  and  (48)  and  stores  the 


computed  wave  amplitudes  Into  the  WAMP  array. 


Subroutines  Called: 


SUBROUTINE  HUXPD 


Function: 

This  routine  computes  the  cosine  and  sine  components  of  wave 
amplitudes  generated  by  point  doublets. 

Calling  Format: 

CALL  HWXPD 

Detailed  Description: 

Subroutine  HWXPD  evaluates  Equations  (45)  and  (46)  for  every  point 
doublet  and  stores  the  computed  wave  amplitude  Into  the  WAMP  array. 

Subroutine  Called: 

None 
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SUBROUTINE  HWXPS 
Function: 

This  routine  computes  the  cosine  and  sine  components  of  wave 
amplitudes  generated  by  point  source-sinks. 

Calling  Format: 

CALL  HWXPS 

Detailed  Descriptions: 

Subroutine  HWXPS  evaluates  Equations  (43)  and  (44)  for  every  point 
source-sink  and  stores  the  computed  wave  amplitudes  Into  the  WAMP  array. 

Subroutines  Called: 

None 
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SUBROUTINE  HUXSD 
Function: 

This  routine  computes  the  cosine  and  sine  components  of  wave 
amplitudes  generated  by  the  surface  doublet.. 

Calling  Format: 

CALL  HWXSD 

Detailed  Description: 

Subroutine  HWXSD  computes  Equations  (39)  and  (41)  for  both  fore  and 
aft  sections  and  stores  the  computed  wave  amplitudes  in  the  WAMP  array. 
The  first  six  rows  of  the  TSC  array  for  the  optimization  of  wavemaking 
resistance  are  also  filled  if  the  value  in  KLP  is  1. 

Subroutine  Called: 


HWXXX 


SUBROUTINE  HUXSS 


Function: 

This  subroutine  computes  the  cosine  and  sine  components  of  the  wave 
amplitudes  generated  by  the  surface  source-sink. 

Calling  Format: 

CALL  HWXSS 

Detailed  Description: 

Subroutine  HWXSS  computes  Equations  (37)  and  (38)  for  each  n-surface 
segment  and  stores  the  wave  amplitudes  computed  in  the  WAMP  array.  The 
first  six  rows  of  the  TSC  array  for  the  optimization  of  wavemaking 
resistance  are  filled  if  the  value  In  KLP  is  1. 

Subroutine  Called: 

HWXXX 
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SUBROUTINE  HWXXX 


Function: 

This  routine  computes  Equations  (35)  and  (36)  or  (40)  and  (42)  for  one 
n-surface  segment. 


Calling  Format:  +  +  +  +  + 

CALL  HWXXX  (XS,  XE,  K,  IND,  L) 


Description  of  Parameters: 


XS.XE  The  starting  and  ending  points  for  integration  interval 
in  the  ^-direction 


The  number  of  n-surface  subsegments 


IND  An  indicator: 

=  1  Surface  source-sink 
=  2  Surface  doublet 


L  An  indicator: 

=  1  The  forebody 
*  2  The  aftbody 


Detailed  Description: 

Subroutine  HWXXX  first  fills  the  arrays  TXP,  R  and  S.  It  then  enters 
into  a  loop  to  compute  Equations  (35)  and  (36)  or  (40)  and  (42),  depending 
on  the  value  of  IND. 


Subroutine  Called: 
PPP 
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SUBROUTINE  HWZSD 


Function; 

T'/;is  r  ur/irki  computes  Equation  (30). 

Cal  inq  Format: 

CALL  HWZSD 

Detailed  Description: 

Subroutine  HWZSD  computes  Equation  (30)  for  both  fore  and  aft 
sections  and  stores  the  results  in  tin  I SD  array. 

Subroutines  Called; 


None 
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PROGRAM  PCWMR 

Function: 

This  program  is  the  main  program  in  the  (2,0)  overlay. 

Calling  Format: 

CALL  AET$KC(5LPCWMR) 

Detailed  Description: 

Program  PCWMR  calls  in  either  HWWMR  or  HWOPT  depending  upon  the 
value  of  KWR.  It  then  calls  in  the  (1,0)  overlay  to  process  other 
input  data. 

Subroutines  Called: 

HWOPT,  HWWMR,  PCIPT 
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FUNCTION  PPP 


Function; 

This  routine  computes  Equations  (35),  (36),  (40)  and  (42)  for  one  i 
value  of  one  n-surface  subsegment. 

Calling  Format:  *  t  t  t  *  ,  *  *****  *  ,  * 

PPP  (CSG.SNG.HAP.HAN,  R1 ,R2,R3,R4,S1 ,S2,S3,S4,TXP1 ,TXP2,TXP3, 
TXP4,A,B,DCX,DCT,IND) 

Description  of  Parameters: 

See  the  program  listing  of  HWXXX 

Detailed  Description: 

Function  PPP  has  two  entry  points,  PPP  and  QQQ;  the  former  for  the 
computation  of  cosine  components,  Equations  (35)  and  (36),  and  the  latter 
for  the  computation  of  the  sine  components.  Equations  (40)  and  (42).  Upon 
entry,  the  value  of  IND  is  checked:  If  IND  =  1,  the  equation  for  the 
surface  source-sink  is  computed;  if  IND  =  2,  the  equation  for  the  surface 
doublet  is  computed. 

Subroutines  Called: 

None 
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SUBROUTINE  HWFCN 
Function: 

This  subroutine  computes  a  coefficient  of  the  wavemaking  resistance. 

Calling  Format:  +  *  ++++++ 

CALL  HWFCN  (PAR,  PART,  N2RET,  PD,  PR,  CM,  CEF,  CR) 

Description  of  Parameters: 

PAR ( 6 )  Six  parameters  for  modifying  the  singularity  distribution 

PART ( 6 )  Six  parameters  for  modifying  the  singularity  distribution 

of  the  last  trial 

N2RET  Indicator  of  search  status: 

=  1  First 

=  2  Subsequent  trial 
=  3  End  of  search  (laKt  trial) 

PD, PR  PD=a,  PR=b  (See  Figure  3) 

CM(6,6)  Refer  to  description  of  CMS  or  CMH  In  Appendix  A 

CEF(6)  Six  coefficients  corresponding  to  the  six  parameters  of 

PAR 

CR  Contains  a  value  computed  by  Eq.  (55) 

Detailed  Description: 

Subroutine  HWFCN  first  converts  the  six  parameters,  PAR(6),  to 
six  coefficients  CEF(6)  (Equation  (25))  and  then  computes  the  coefficient 
(CR)  of  wavemaking  resistance. 

Subroutine  Called: 

HWPTC 
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SUBROUTINE  HWOLD 


Function: 

This  routine  determines  the  six  parameters  for  the  horizontal  line 
doublet  of  either  forebody  or  aftbody  that  produce  the  lowest  wavemaking 
resistance. 

Calling  Format: 

CALL  HWOLD 

Detailed  Description: 

Subroutine  HWOLD  makes  successive  calls  to  Subroutine  EXPRS  to  randomly 
determine  a  set  of  six  parameters  corresponding  to  the  horizontal  line 
doublet  (Equation  (58)).  In  each  successive  call.  Subroutine  HWFCN  com¬ 
putes  the  wavemaking  resistance  using  the  six  parameters  generated. 
Subroutine  EXPRS  then  compares  the  wavemaking  resistance  value  just  com¬ 
puted  with  the  lowest  value  obtained  thus  far.  If  the  newly  computed 
value  is  lower,  the  six  parameters  used  to  produce  this  lower  value  will 
replace  the  six  parameters  previously  saved.  At  the  end  of  the  search, 
the  six  parameters  saved  (those  parameters  which  produced  the  lowest 
wavemaking  resistance  value)  will  be  converted  Into  their  coefficient 
values  and  these  values  will  be  added  to  the  values  of  the  original 
horizontal  line  doublet. 

Subroutines  called: 

EXPRS,  HWCTP ,  HWFCN,  HWPTS 
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SUBROUTINE  HWOLS 


Function: 

This  routine  determines  the  six  parameters  for  the  horizontal 
line  source-sink  of  either  forebody  or  aftbody  that  produce  the  lowest 
wavemaking  resistance. 

V 

Calling  Format: 

CALL  HWOLS 

Detailed  Description: 

Subroutine  HWOLS  makes  successive  calls  to  Subroutine  EXPRS  to 
randomly  determine  a  set  of  six  parameters  corresponding  to  the  horizontal 
line  source-sink  (Equation  (58)).  In  each  successive  call,  Subroutine 
HWFCN  computes  the  wavemaking  resistance  using  the  six  parameters 
generated.  Subroutine  EXPRS  then  compares  the  wavemaking  resistance 
value  just  computed  with  the  lowest  value  obtained  thus  far.  If  the 
newly  computed  value  is  lower,  the  six  parameters  used  to  produce  this 
lower  value  will  replace  the  six  parameters  previously  saved.  At  the 
end  of  the  search,  the  six  parameters  saved  (those  parameters  which 
produced  the  lowest  wavemaking  resistance  value)  will  be  converted  into 
their  coefficient  values  and  these  values  will  be  added  to  the  values 
of  the  original  horizontal  line  source-sink. 


Subroutines  called: 

EXPRS,  HWCTP ,  HWFCN,  HWPTS 
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SUBROUTINE  HWOPT 
Function: 

This  subroutine  is  the  executive  program  for  the  minimization  of 
wavemaking  resistance. 

Calling  Format: 

CALL  HWOPT 

Detailed  Description: 

Subroutine  HWOPT  reads  in  the  data  for  the  optimization  process. 
Using  the  data  provided,  the  program  optimizes  one  singularity  distribu¬ 
tion  for  either  the  forbody  or  the  aftbody  using  a  single  J  value. 

It  then  calls  HWWMR  to  store  data  into  the  array  TSC,  after  which  it 
calls  one  of  the  following  subroutines  to  evaluate  the  six  coefficients 
of  Equation  (58):  HWOSS  (for  surface  source-sink),  HWOSD  (for  surface 
doublet),  HWOLS  (for  horizontal  line  source-sink),  or  HWOLD  (for 
horizontal  line  doublet).  For  detailed  description,  refer  to  Figure  4. 

Subroutine  called: 

HWOLS,  HWOLD,  HWOSS,  HWOSD,  HWWMR 
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SUBROUTINE  HWOSD 


Function: 

This  routine  determines  the  six  parameters  for  the  surface  doublet 
of  either  forebody  or  aftbody  that  produce  the  lowest  wavemaking 
resistance. 

Calling  Format: 

CALL  HWOSD 

Detailed  Description: 

Subroutine  HWOSD  makes  successive  calls  to  Subroutine  EXPRS  to 
randomly  determine  a  set  of  six  parameters  corresponding  to  the  surface 
doublet  (Equation  (58)).  In  each  successive  call.  Subroutine  HWFCN 
computes  the  wavemaking  resistance  using  the  six  parameters  generated. 
Subroutine  EXPRS  then  compares  the  wavemaking  resistance  value  just 
computed  with  the  lowest  value  obtained  thus  far.  If  the  newly  computed 
value  is  lower,  the  six  parameters  used  to  produce  this  lower  value 
will  replace  the  six  parameters  previously  saved.  At  the  end  of  the 
search, the  six  parameters  saved  (those  parameters  which  produced  the  lowest 
wavemaking  resistance  value)  will  be  converted  into  their  coefficient 
values  and  these  values  will  be  added  to  the  values  of  the  original 
surface  doublet. 


Subroutines  called: 


EXPRS,  HWCTP,  HWFCN,  HWPTS 


SUBROUTINE  HWOSS 


Function; 

This  routine  determines  the  six  parameters  for  the  surface  source- 
sink  of  either  forebody  or  aftbody  that  produce  the  lowest  wavemaking 
resistance. 

Calling  Format; 

CALL  HWOSS 

Detailed  Description; 

Subroutine  HWOSS  makes  successive  calls  to  Subroutine  EXPRS  to 
randomly  determine  a  set  of  six  parameters  corresponding  to  the  surface 
source-sink  (Equation  (58)).  In  each  successive  call.  Subroutine  HWFCN 
computes  the  wavemaking  resistance  using  the  six  parameters  generated. 
Subroutine  EXPRS  then  compares  the  wavemaking  resistance  value  just 
computed  with  the  lowest  value  obtained  thus  far.  If  the  newly  computed 
value  Is  lower,  the  six  parameters  used  to  produce  this  lower  value  will 
replace  the  six  parameters  previously  saved.  At  the  end  of  the  search, 
the  six  parameters  saved  (those  parameters  which  produced  the  lowest 
wavemaking  resistance  value)  will  be  converted  Into  their  coefficient 
values  and  these  values  will  be  added  to  the  values  of  the  original 
surface  source-sink. 

Subroutines  called; 

EXPRS,  HWCTP,  HWFCN,  HWPTS 
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SUBROUTINE  OAUX 
Function: 

This  routine  directs  the  calling  of  other  routines  to  calculate 
the  u,  v,  w  contributions  from  the  various  singularity  distributions. 

Calling  Format:  +  + 

CALL  DAUX  (A,  B,  F) 

Description  of  Parameters: 

A  Contains  the  x-value  of  a  streamline  point  being  evaluated 

B(2)  Contains  the  y-  and  z-values  of  a  streamline  point  being 

evaluated 

F(2)  Contains  slopes  of  the  streamline  points  being  evaluated 
Detailed  Description: 

Subroutine  DAUX  checks  on  each  singularity  distribution  calculation 
flag  periodically  to  determine  If  the  contribution  from  that  particular 
type  of  distribution  Is  to  be  calculated.  If  It  is,  a  call  Is  made  to 
the  appropriate  routine.  When  the  calculation  is  completed,  4n  is  added 
to  the  U  component  (that  is,  the  normalized  steady  flow  V  =  4u)  and  the 
slopes  of  the  streamlines  are  calculated  from  the  streamline  formula. 
Equation  (10). 

Subroutines  Called: 


SURF,  PSOUR,  PDOUB,  LSOURV,  LDOUBV,  LSOURH,  LDOUBH 


SUBROUTINE  INT 


Function: 

This  routine  performs  an  Integration  of  one  variable  of  a 
specific  family  of  functions  over  some  definite  limits. 


Calling  Format: 


t  + 


CALL  INT  (XXI,  YET,  Z,  F,  G,  S,  D) 


Description  of  Parameters: 
F 


s(,)  ’/ 


X1'1  dx 


rG  (X2-  2-Z-X  +  XXI2  +  YET2  +  z2)3/2 


I  *  1,2,. ...7 


F  X1"1  dx 

D(I)  sf  (X^-  2-Z-X  +  XX 1 2  +  YET2  ♦  Z2)5/Z  1  *  1  *2 . 8 

G 

Detailed  Description: 

In  order  to  save  multiplications  and  other  arithmetic  operations 
later,  subroutine  INT  begins  by  calculating  several  frequently  used 
Intermediate  quantities.  Then,  by  checking  the  various  distribution 
calculation  flags,  those  values  (and  only  those  values)  of  S  and  D 
needed  by  the  calling  routine  are  determined  by  using  the  limits  in  the 
analytical  solutions  of  the  above  integrals. 


Subroutines  Called: 
None 
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SUBROUTINE  LDOUBH 


Function: 

This  routine  calculates  the  u,  v,  w  contribution  from  all  horizontal 
line  doublets. 

Calling  Format:  +  +  + 

CALL  LDOUBH  (ULDH,  VLDH,  WLDH) 

Description  of  Parameters: 

ULDH  I 

VLDH  I  u,  v,  w  components,  respectively,  of  the  streamline 
WLDH  ,elocU* 

Detailed  Description: 

Subroutine  LDOUBH  performs  the  same  operations  on  horizontal -line 
doublets  that  LSOURH  performs  on  horizontal -line  source-sinks  (See 
Equations  (95)  through  (97).) 

Subroutines  Called: 

I  NT 


SUBROUTINE  LDOUBV 


Function: 

Tnis  routine  calculates  the  u,  v,  w  contributions  from  all 
vertical -line  doublets. 


Calling  Format:  +  + 

CALL  LDOUBV  (ULOV,  VLDV,  WLDV) 


Description  of  Parameters: 


ULDV 

VLDV 

WLDV 


u,  v,  w  components,  respectively,  of  the  streamline 
velocity 


Detailed  Description: 

Subroutine  LDOUBV  performs  the  same  operations  on  vertical-line 
doublets  that  LSOURV  performs  on  vertical-line  source-sinks  (See 
Equations  (87)  through  (89).) 


) 


I 


i 

i 

! 
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SUBROUTINE  LSOURH 


Function: 

This  routine  calculates  the  u,  v,  w  contributions  from  all 
horizontal  line  source-sinks. 


Detailed  Description: 

After  setting  the  distribution  calculation  flag  to  active  (l.e., 
making  KLSH  <0),  subroutine  LSOURH  calculates  YY  and  11.  The  Integration 
over  the  line  source-sink  is  performed  and  the  results  used  to  calculate 
the  contributions.  Each  line  Is  also  checked  to  see  if  it  is  off  the 


centerline  and/or  below  the  water  plane.  If  it  Is,  the  corresponding 
appropriate  Image  line  source-sinks  are  deactivated  by  being  made 
positive  again.  (See  Equations  (91)  through  (93).) 


Subroutines  Called: 


SUBROUTINE  LSOURV 


Function: 

This  routine  calculates  the  u,  v,  w  contributions  from  all  vertical 
line  source-sinks. 


Calling  Format: 


+ 


CALL  LSOURV  (ULSV,  VLSV,  WLSV) 


Description  of  Parameters; 


u,  v,  w  components,  respectively,  of  the  streamline 
velocity. 


Detailed  Description: 

After  setting  the  distribution  calculation  flag  to  active  (by 
making  KLSV<0),  subroutine  LSOURV  calculates  XX,  YY  and  11.  The 
Integrations  over  the  line  source-sink  and  Its  Image  are  performed  and 
the  results  used  to  calculate  the  contributions.  When  the  line  Is  not  on 
the  centerline,  an  Image  source-sink  distribution  Is  created  on  the 
other  side  and  the  contribution  from  that  Image  and  Its  image  above  the 
waterplane  is  calculated.  KLSV  Is  then  deactivated  by  being  made 
positive  again.  (See  Equations  (83)  through  (85).) 


Subroutines  Called: 
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SUBROUTINE  PCSTL 


Function: 


This  routine  Is  the  main  program  In  the  (3,0)  overlay. 


Calling  Format: 

CALL  AETSKC  (5LPCSTL) 


Detailed  Description: 

After  Subroutine  PCSTL  Initializes  several  variables.  It  calls 
STREAM  to.  trace  the  streamlines  and  then  calls  in  the  (1,0)  overlay 
to  process  other  input  data. 


Subroutines  Called: 


STREAM,  PCIPT 
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SUBROUTINE  PDOUB 


Function: 

This  routine  calculates  the  u,  v,  w  contribution  from  all  point 
doublets. 

Calling  Format:  +  +  + 

CALL  PDOUB  (UPOD,  VPDD,  WPDD) 


Description  of  Parameters: 


UPDD 

VPDD 

WPDD 


u,  v,  w  components,  respectively,  of  the  streamline 
velocity. 


Detailed  Description: 

Subroutine  PDOUB  performs  the  same  operations  on  point  doublets 
that  PSOUR  performs  on  point  source-sinks.  (See  Equations  (79)  through 


t 

(81).) 


Subroutines  Called: 
None 
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SUBROUTINE  PSOUR 


Function; 

This  routine  calculates  the  u,  v,  w  contributions  from  all  point 
source-sinks. 


Calling  Format;  +  +  + 

CALL  PSOUR  (UPSD,  VPSD,  WPSD) 


Description  of  Parameters; 

UPSD 

VPSD  .  u,  v,  w  components,  respectively,  of  the  streamline  velocity 
WPSD 


Detailed  Description; 

Each  point  source-sink  is  examined  to  see  if  it  is  off  the  center- 
line  and/or  below  the  waterplane.  If  it  is,  appropriate  image  source- 
sinks  are  created.  The  contributions  from  the  source-sink  and  any 
images  are  then  calculated.  (See  Equations  (75)  through  (77).) 


Subroutines  Called; 


None 


SUBROUTINE  QUAD 


Function: 

Using  a  10-point  Gaussian  quadrature  method,  this  routine  numerically 
evaluates  the  double  Integrals  used  In  calculating  the  u,  v,  w  contri¬ 
butions  from  the  surface  distributions. 

Calling  Format:  +  +  +  + 

CALL  QUAD  (XL,  XU,  S,  IFL) 

Description  of  Parameters: 

XL, XU  Limits  of  Gaussian  Interval  being  evaluated  over 

S (3)  Array  containing  the  u,  v,  w  contributions 

IFL  Flag  indicating  which  side  of  n-surface  Is  being  evaluated. 

IFL>0,  n  >  0;  IFL>  0,  n  <  0. 

Detailed  Description: 

This  routine  evaluates  the  six  double  integrals  In  Equations  (67), 

(68),  (69),  (71),  (72)  and  (73). 

The  entire  integration  is  performed  so  as  to  avoid  any  unnecessary 
duplication  of  arithmetic  computation  and  thus  to  save  time.  The  actual 
calculation  of  the  functions  and  the  summing  of  the  results  for  these  ten 
points  is  performed  separately  for  the  source  and  doublet  cases,  with  the 
results  being  stored  in  the  array  D.  D(l)  contains  the  u  component  of  the 
velocity,  D(2)  the  v,  and  D(3)  the  w.  These  results  are  then  multiplied 
by  the  appropriate  distribution  coefficient  and  summed  into  the  TT  array. 

Finally,  a  last  postmultiplication  is  performed  on  TT  corresponding  to  the 
multiplication  by  b-a  In  the  quadrature  formula.  The  final  results  are 
summed  Into  S  for  return  to  SURF. 

Subroutines  Called: 

I  NT 

1 
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SUBROUTINE  STREAM 


Function: 

This  routine  acts  as  the  executive  for  the  actual  streamline  calculations. 

Calling  Format: 

CALL  STREAM 

Detailed  Description: 

Subroutine  STREAM  first  reads  In  the  number  of  stations,  the  number  of 
streamlines,  the  station-number  option,  the  offset  scaling  factor,  and 
other  variables.  It  then  reads  the  x-value  of  streamline  stations  and 
the  Initial  y  and  z  offsets  for  each  streamline.  The  Kutta-Merson 
routine  KUTMER  Is  then  called,  a  loop  being  provided  each  time  with  the 
x,  y,  z  coordinates  of  the  present  streamline  point  and  the  a  x  to  the 
next  station.  The  routine  KUTHER  returns  the  x,  y,  z  coordinates  of  the 
next  point.  These  points  are  saved  and  then  printed  out  when  the  last 
station  Is  reached  or  when  KUTMER  falls  to  meet  the  error  criteria 
specified  In  EPS  and  A. 

Subroutines  Called: 

KUTMER 
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SUBROUTINE  SURF 


Function: 

This  routine  calculates  the  u,  v,  w  contribution  from  the  surface 
source-sink  and  doublet  distributions. 


Calling  Format:  +  +  + 
CALL  SURF  (A,  B,  C) 


Description  of  Parameters: 

A 

B  u,  v,  w  components,  respectively,  of  the  streamline  velocity. 

C 

Detailed  Description: 

Subroutine  SURF  Initializes  IFL  and  then  calculates  the  value  of  ET. 

Y  Is  checked  to  see  If  It  has  penetrated  the  n-surface.  In  which  case  a 

message  Is  printed  out.  XINT  Is  then  calculated  and  preparations  are  made 

to  begin  the  Integration  on  the  +y  side  of  the  n-surface  In  the  +x  direction. 

XL  and  XU  are  calculated,  with  XL  being  the  last  value  of  XX  and  XU  being 
2 

XX  +  XINT  *  J  .  Starting  at  1,  0  Increases  by  1  for  each  Interval.  When 
J  >  5,  the  remaining  portion  of  the  n-surface  Integration  In  that  direction 
Is  computed  In  one  final  step.  The  values  of  XL  and  XU  are  checked  to  see 
that  they  (1)  stay  within  the  n-surface  limits,  (2)  lie  outside  of  the 
parallel  body  lengths,  and  {3}  both  lie  on  the  same  fore  and  aft  section. 
Changes  are  made  to  XL,  XU  and/or  XX  where  appropriate  to  meet  the  above 
conditions.  Then  If  XL  >  XU,  a  call  Is  made  to  the  Gaussian  quadrature 
routine  QUAD  which  will  calculate  the  actual  contribution  ovet  the 
Interval  XL+XU.  This  process  Is  repeated  until  the  Integration  Interval 
reaches  the  forward  end  of  the  n-surface,  x  =  1.0.  The  whole  procedure  Is 
then  repeated  to  perform  the  Integration  from  X  In  the  -x  direction 
and  then  again  to  do  the  Integration  In  both  directions  on  the  -n  side  of 
the  n-surface. 

Subroutines  Called: 

QUAD 


SUBROUTINE  VOLAR 


Function: 

This  routine  computes  the  value  and  wetted  surface  for  the  strut 
or  lower  hull  for  SWATH  ship. 

Calling  Format: 

CALL  VOLAR  (OFF,  ZS,  KPT,  YH,  ZZZ,  KXS,  KXE) 

Description  of  Parameters: 

OFF  Offset  scaling 

ZS  Depth  of  the  strut  or  lower  hull 

KPT  Station  number  options: 

»1  Traces  streamline  at  every  station 
®2  Traces  streamline  at  every  other  station 
YM  Beam  of  the  strut  or  the  lower  hull 

ZZZ  Draft  of  the  strut 

KXS  Station  number  that  marks  the  fore  point  of  the  bow 
at  the  waterline 

KXE  Station  number  that  marks  the  aft  point  of  the  stern 
at  the  waterline 

Detailed  Description: 

Subroutine  VOLAR  calls  subroutine  SIMPUN  to  compute  the  volume  and 
wetted  surface  for  the  strut  or  lower  hull  of  the  SWATH  ship. 

Subroutine  called: 


SIMPUN 
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APPENDIX  A 
VARIABLES  IN  COMMON 

The  variable  names  used  In  COMMON  throughout  the  program  are  described 
In  the  pages  that  follow.  The  User's  Manual  contains  detailed  descriptions 
of  the  Input  data  variables.  If  the  Input  data  variables  are  also 
Included  In  COMMON,  only  the  names  of  the  Input  data  variables  and  the 
corresponding  numbers  of  the  cards  containing  their  descriptions  are 
Indicated. 


AAA(N,L,K)  The  value  of  A^  used  In  Equation  (13)  to  approximate  the 
n-surface 

N  The  number  of  straight-line  segments  (N<4)  defining  the 
n-surface 


L 


K 


AKZ 

BBB(N,L,K) 


BEM 

BEXT(N) 

BKZ 


The  portion  of  the  ship  under  consideration: 

Forebody 
=2  Aftbody 

The  singularity  distribution  to  be  calculated: 

*1  Surface  source-sink 
=2  Surface  doublet 

=3  Reserved  for  Interactive  graphics  use 
=4  Reserved  for  interactive  graphics  use 

A  constant  defined  as  j^j-  where  KQ  =  9/y2 

The  value  of  A2  used  In  Equation  (13)  to  approximate  the 
n-surface.  The  values  for  N,L,  and  K  are  the  same  as 
those  for  AAA(N,L,K)  above 

[Defined  on  Input  data  card  C8] 

Values  of  exponents  for  exponential  random  search,  where  N=3 
A  constant  defined  a  KQ  tane  sece,  where  K0=9/v2 


f 


Preceding  pigi  blink 
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CESD(I,J,L) 

L 


CESS(I,J,L) 

L 


CLDH(I,L) 

I 


CLDH(7, 


Coefficients  of  Equation  (16)  for  the  surface  doublet, 
where  1  <^I<j6  and  1  <J<4 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

=3  Reserved  for  Interactive  graphics  use 
=4  Reserved  for  interactive  graphics  use 


Coefficients  of  Equation  (16)  for  the  surface  source-sink, 
where  1  <_I<_6  and  1  <J<4 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

=3  Reserved  for  Interactive  graphics  use 
=4  Reserved  for  interactive  graphics  use 

Information  defining  horizontal -line  doublet 


Indicator: 

!  Coefficients  of  Equation  (20) 

7  j  Integration  interval  in  the  ^-direction 

A 


2) 


CLDH(8,2) 


CLDH(7,1 ) 


CLDH(8,1 ) 


=9  The  distance  from  the  centerline  of  the  line  doublet 
to  the  ?-£  plane 

=10  The  distance  from  the  center  of  the  line  doublet  to 
the  s-n  plane 


L 


Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

=3  Reserved  for  interactive  graphics  use 
=4  Reserved  for  interactive  graphics  use 

CLDV(J,L)  Information  defining  vertical  line  doublet 

J  Indicator: 

=1  \ 

•  V 

./Coefficients  of  Equation  (19) 

4) 

=5  Contains  distance  from  the  centerline  of  the 
line-doublet  to  the  n-c  plane 
=6  Contains  distance  from  the  centerline  of  the 
line-doublet  to  the  plane 

=8 }  Contains  integration  interval  in  the  ^-direction 


L  Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

=3  Reserved  for  interactive  graphics  use 
=4  Reserved  for  interactive  graphics  use 
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CLSH(I,L) 

I 


L 


CLSV(J,L) 

J 


i 

j 


"T7T>"  «T  -rfrty 


Information  defining  horizontal  line  source-sink 
Indicator: 

!  Coefficients  of  Equation  (20) 

~g|  Integration  interval  in  the  ^-direction 


A' 


* 

CLSH(7,2)  CLSH(8,2) 

:LSH(7,1)  CLSH(8,1) 

=9  Distance  from  centerline  of  the  line  source-sink 
to  the  s-5  plane 

=10  Distance  from  center  of  the  line  source-sink  to 
the  5-n  plane 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

=3  Reserved  for  interactive  graphics  use 
=4  Reserved  for  interactive  graphics  use 

Information  defining  the  vertical  line  source-sink 

Indicator: 

i  Contains  coefficients  for  Equation  (19) 

=5  Distance  from  the  centerline  of  the  line  source- 
sink  to  the  n-c  plane 

=6  Distance  from  the  centerline  of  line  source-sink 
to  the  ?-£  plane 

“gj  Integration  interval  in  the  c-direction 


CLSV(8,2) 

CLSV(7,2) 


CLSV (8,1 ) 
CLSV (7,1 ) 


L 


CMH(M,N,L) 

L 


CMS(M,N,L) 

L 


CN3(N) 

CRW(N) 

CSN 

DDR 

DEG(N) 


Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

=3  Reserved  for  Interactive  graphics  use 
=4  Reserved  for  Interactive  graphics  use 

Coefficients  of  the  Matrix  (24)  for  horizontal -line 
source-sink  and  doublet,  where  M=7  and  N=7 

Hull  portion  to  be  considered: 

=1  Forebody,  horizontal  line  source-sink 
=2  Aftbody,  horizontal  line  source-sink 
=3  Forebody,  hor  .ontal  line  doublet 
=4  Aftbody,  horizontal  line  doublet 

Coefficients  of  the  Matrix  (24)  for  surface  source-sink 
and  doublet,  where  M=7  and  N=7 

Hull  portion  to  be  considered: 

=1  Forebody,  surface  source-sink 
=2  Aftbody,  surface  source-sink 
=3  Forebody,  surface  doublet 
=4  Aftbody,  surface  doublet 

3 

Constants  defined  as  Cos  0  where  1<N<NDG 

Wavemaking  resistance  Equation  (55)  at  different  Froude 
numbers,  where  1<N<NDG 

A  constant  defined  as  Cose 

[Defined  on  input  data  card  C24] 

The  values  in  degree  of  direction  angle  where  1<N<NDG  and 
0<DEG(N)<85° 
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OGR ( N )  Information  for  numerical  integration  of  Equation  (55) 

N  Indicator: 

=1  Initial  increment  size,  in  degrees 

=2  Upper  limit  of  the  integration  where  0<D6R(2)<85° 

=3) 

•/Those  points  (degrees)  at  which  the  increment  size 
‘  (is  to  be  reduced  by  half 
1CK 


DIN  The  increment  size  of  0  for  current  use  in  the  numerical 

integration  of  Equation  (55) 

DXX(L,N)  Integration  subinterval  for  surface  source-sink  doublet 


FOE 

[Defined  on  input  data  card  C3] 

FDI 

[Defined  on  input  data  card  C3] 

FDP(N) 

[Defined  on  input  data  card  C75] 

FDS 

[Defined  on  input  data  card  C3] 

FDW(N) 

[Defined  on  input  data  card  C4] 

FND 

Variable  NDY  in  floating-point  format 

FRD(N) 

Froude  numbers  at  which  wavemaking  resistance  will  be 
evaluated,  1<N<MFD 

KAYT(N) 

[Defined  on  input  data  card  C72] 

KBD 

[Defined  on  input  data  card  C5] 

KBS 

[Defined  on  input  data  card  C5] 

KED(J.L) 

Indicators  for  the  calculation  of  coefficients  of  the 

surface  doublet 
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J 

An  Indicator: 

=0  CESD(I,J,L)  =  0,  1=  1 .6 
=1  CESD(I,J,L)  t  0,  1=1,. ..6 

L 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

KES(J,L) 

Indicators  for  the  calculation  of  coefficients  of  the 

surface  source-sink: 

J 

An  indicator: 

=0  C£SS(I ,J,L)  =  0,  I  =  1 , . . .6 
*1  CESS(I ,J,L)  t  0,  I  =  1 ... .6 

L 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

KLDH 

[Defined  on  input  data  card  C5] 

KLDV 

[Defined  on  input  data  card  C5] 

KLP 

An  indicator  used  in  the  optimization  process: 

=1  Forebody 
=2  Aftbody 

KLSH 

[Defined  on  input  data  card  C5] 

KLSV 

[Defined  on  input  data  card  C5] 

KOP(K) 

[Dpfined  on  input  data  card  C73] 

KPD 

[Defined  on  input  data  card  C5] 

KPI(K) 

[Defined  on  input  data  card  C7] 

KPR(K) 

[Defined  on  input  data  card  C6] 

KPS 

[Defined  on  input  data  card  C5] 

KSD 

[Defined  on  input  data  card  C5] 

KSS 

[Defined  on  input  data  card  C5] 

KWH  An  Indicator  used  In  the  optimization  process: 

s0  No  optimization 
*1  Optimization  for  the  forebody 
*-l  Optimization  for  the  aftbody 

KWR  [Defined  on  input  data  card  C2] 

KWZ  The  current  KZP  value 

KZP(N)  [Defined  on  input  data  card  C74] 

LING  An  Indicator: 

=0  Standard  option  in  which  the  first  random  number 

chosen  is  compared  to  .5  to  determine  the  search 
region  ft  r  each  trial 

=1  The  random  number  is  compared  with  a  variable  Instead  of 
with  a  constant  of  .5 

MDG  A  constant  which  is  the  maximum  dimension  for  the  arrays 

DEG  and  THE,  where  MDG  *  100 

MEND  The  number  of  stations  in  the  streamline  tracing 

MFD  A  constant  which  Is  the  maximum  number  of  Froude  numbers 

allowed,  where  MFD  =  30 

NDG  The  index  to  the  arrays  DEG  and  THE;  1<NDG£MDG 

NDY  [Defined  on  Input  data  card  C8] 

NFD  The  index  to  the  array  FRD.  1<NFD<MFD 

NFP  [Defined  on  input  data  card  C75] 

NFW  The  index  to  the  array  FDW.  1<NFW<MFW 

NGR  A  constant  which  is  the  number  of  times  the  variable  DIN 

will  be  divided  by  2 

NIT  [Defined  on  input  data  card  C8] 

NPD  The  total  number  of  point-doublets  read  from  the  input 

file 
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NPS 


The  total  number  of  point-source-sinks  read  from  the  input 
file 


NTI 

A  constant  which  is  the  i  va  ue  in  the  Equation  (15), 
where  NTI=6 

NTJ 

A  constant  which  is  the  j  valu?  in  the  Equation  (15), 
where  NTJ=4 

NTS 

A  constant,  hwere  f.TS=8 

NZP 

[Defined  on  input  data  card  C74] 

N2EXP 

Number  of  exponents,  where  N2EXP=3 

N2PAR 

Number  of  search  parameters,  where  N2PAR=6 

PALS(I.N.L) 

Parameters  of  the  singularity  distribution  of  the 
horizontal  line  source-sink  and  doublet 

I 

Indicator: 

.  >The  six  parameters 

N  An  indicator: 

=1  Used  in  the  interactive  graphics 
=2  Currently  used 

L  Hull  portion  to  be  considered: 

=1  Forebody,  horizontal  line  source-sink 
=2  Aftbody,  horizontal  line  source-sink 
=3  Forebody,  horizontal  line  doublet 
=4  Aftbody,  horizontal  line  doublet 

PARD(I,N,L)  Parameters  of  the  singularity  distribution  of  the  surface 

doublet 

I  Indicator: 

=1) 

:  (The  six  parameters 
6 
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N 

L 

PARH(6) 

PARL(G) 

PARS(I,N,L) 

I 


N 

L 

PCY 

PIE 

PIND(K,N) 
PINS (N ,N ) 
PMD(L) 


Indicator: 

=1  Reserved  for  Interactive  graphics 
-2  Currently  used 

Hull  portions  to  be  considered: 

=1  Forebody 
-2  Aftbody 

[Defined  on  Input  data  card  C79] 

[Defined  on  Input  data  card  C78] 

Parameters  for  the  singularity  distribution  of  the  surface 
source-sink 

Indicator: 

The  six  parameters 

Indicator: 

=1  Reserved  for  interactive  graphics 

-2  Currently  used 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

[Defined  on  input  data  card  C2] 

A  constant  which  is  equal  to  n 
[Defined  on  input  data  card  C42] 

[Defined  on  input  data  card  C40] 

The  lengths  of  the  paralle’  mid-body  for  surface  doublet 
computation 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 
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PMS(L) 

L 


PS(L) 

L 


RAD 

RKZ 

RLM 

SNE 

T 

TEF(I.J) 

THE(N) 

TSC(M,N,K) 


WAMP(N,L,K) 

L 


The  lengths  of  the  parallel  midbody  for  surface  source- 

sink  computation 

Hull  portion  to  be  considered 

=1  Forebody 

=2  Aftbody 

The  lengths  of  the  parallel  midbody  of  either  surface 
source-sink  or  doublet 

Hull  portion  to  be  considered: 

=1  Forebody 
=1  Aftbody 

The  constant  57.2958 

The  variable  which  is  defined  ko  where  ko  =  9/y2 
The  variable  KQ.SEC^e 
The  variable  Sine 
Same  as  for  ZRH 

A  temporary  array  used  tu  store  coefficients  of  singularity 
distributions  in  the  optimization  of  wavemaking  resistance 

Direction  angles  in  radius,  .'here  1<N<100 

Coefficients  of  Matrices  (63)  and  (65) 

Indicator: 

=1  Matrix  (63) 

=2  Matrix  (65) 

The  wave  amplitudes  at  different  direction  angles,  where 
1<N<NDG 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 
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K  Indicator: 

=1  Surface  source-sink 
=2  Surface  doublet 
=3  Point  source-sink 
=4  Point  doublet 
=5  Vertical  line  source-sink 
=6  Vertical  line  doublet 
=7  Horizontal  line  source-sink 
=8  Horizontal  line  doublet 


WTP(N)  [Defined  on  Input  data  card  C77] 

X  X- values  of  a  streamline  point  at  which  evaluation  of  u, 

v,  w  Is  being  made 


XIE(L,K)  Integration  intervals  for  surface  source-sink  and  doublet 

L  Hull  portion  to  be  considered: 

*1  Forebody 
=2  Aftbody 
K  Indicator: 


=1  Surface  source-sink 
=2  Surface  doublet 


i 

: 

< 

XIS(2,*) 

XIE(2,*) 

xis(i>) 

XIE(l,*f 

XIS(L,K)  The  integration  intervals  for  surface  source-sink  and 

doublet  (see  preceding  figures) 

L  Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 
K  Indicator: 

=1  Surface  source-sink 
=2  Surface  doublet 
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Integration  intervals  for  the  surface  doublet 


Indicator: 

=1  Starting  point 

=2  Ending  point 

Hull  portion  to  be  considered: 

=1  Korebody 

=2  Aftbody  _ 


A 


XSD(2,2) 


XSD(1 ,1 ) 


XSD(2,1 ) 


Equation  (35),  where  1=1 ,2 . 6 

Equation  (36),  where  1-1,2..., 6 

Integration  Intervals  for  the  surface  source-sink 
Indicator: 

=1  Start  point 
=2  End  point 

Hull  portion  to  be  considered: 

=1  Forebody 
=2  Aftbody 

t 


>S(1 ,2) 


- 

XSS(2,2) 


- >  - 

XSS(1,1) 


— • — c 

XSS (2,1 ) 


X-values  of  streamlines 

Number  of  streamline  stations  where  1<M<61 

Number  of  streamlines  where  1<N<6 

Y-values  of  a  streamline  point  at  which  evaluation  of  U, 
V,  W  is  being  made 
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YY(M,N)  Y-values  of  streamlines 

M  Number  of  streamline  stations  1<M<61 

N  Number  of  streamlines,  1<N<6 

Z  Z-values  of  a  streamline  point  at  which  evaluation  of 

U,V,W  Is  being  made 

ZER  A  constant  set  to  0.0000001 

ZLS  Result  of  Equation  (30),  for  horizontal  line  source- 

sink  or  doublet 

ZRH  [Defined  on  input  data  card  C8] 

ZSD(J)  Result  of  Equation  (30),  for  surface  source-sink  or 

doublet 

ZZ(M,N)  Z-values  of  streamlines 

N  Number  of  streamline  stations,  1<M<61 

N  umber  of  station-lines,  1<N<6 
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APPENDIX  B 


COMPARISONS  OF  COEFFICIENTS  OF  THEORETICAL  WAVEMAKING 
RESISTANCE  AND  RESIDUAL  RESISTANCE  OF 
SWATH'S  III,  IV,  V,  AND  RC-2 

In  this  appendix,  resistance.,  derived  from  experimental  tests  of 
certain  SWATH  models  are  compared  with  wavemaking  resistances  predicated 
for  the  SWATH  designs  by  the  Plen  program.  Four  recently  developed  SWATH 
models  are  used  In  the  comparison:  SWATH's  III,  IV,  and  V,  developed  at 
the  Naval  Ship  Research  ^nd  Development  Center,  and  SWATH  RC-2,  developed 
at  the  Naval  Undersea  Center.  The  SWATH-model  experimental  tests  were 
Informally  documented  within  the  Center  conducting  the  experiment. 

The  coefficient  of  the  wavemaking  resistance  Cy  Is  defined  as 

y  (1/2)  p  •  S  •  V 

where  Ry  is  the  wavemaking  resistance  In  pounds 
p  is  the  mass  density 

S  is  the  wetted  surface  for  both  strut  and  lower  hull 
V  is  the  ship  speed  In  feet  per  second 

Since  Lg  was  selected  as  the  effective  ship  length  for  the  experimental 
tests,*  Le  was  used  in  the  Pien  program  also  for  the  computation  of  the 
speed-length  ratio  V//L^  where  V  is  the  ship  speed  In  knots.  The 
limitations  of  wavemaking  resistance  theory  are  well  knrwn;  hence  several 


Le* 


where 


!s  •  V  sh  •  'h 

'V  V 

Ls  is  the  length  of  the  strut 

L^  Is  the  length  of  the  lower  hull 

S$  is  the  wetted  surface  of  the  strut 

S^  is  the  wetted  surface  of  the  lower  hull 
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from  a  singularity  distribution  as  discussed  In  Pien's  paper  on  "Motion 

O 

and  Resistance  of  a  Low-Waterplane-Area  Catamaran."  Nevertheless,  as 


the  Figures  10  to  37  which  compare  the  Cr  (residual  resistance  coefficient) 
and  the  Cy  (wavemaking  resistance  coefficient)  show,  the  Plen  program  is 
more  than  adequate  for  estimating  the  resistance  for  a  proposed  SWATH 


design.  This  program  is  also  a  useful  tool  for  making  systematic 


parametric  studies  of  SWATH  ships. 
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Figure  15  -  Ship  and  Experimental  Model  Test  Data  for  SWATH  III  (Model  5276E) 
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Figure  19-  Abbreviated  Lines  and  Profile  for  SWATH  IV  Represented  by  Model  5287 


Figure  20  -  Abbreviated  Body  Plan  for  SWATH  IV  (Model  5287) 
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Figure  21  -  Ship  and  Experimental  Model  Test  Data  for  SWATH  IV  (Model  5287) 


Figure  22  -  Residual  and  Wavemaking  Resistance  Coefficients  versus  Speed-Length  Ratio 
for  SWATH  IV  (Model  5287)  at  a  32-Foot  Draft 


Figure  23  -  Residual  and  Wavemaking  Resistance  Coefficients  versus  Speed-Length  Ratio 
for  SWATH  IV  (Model  5287)  at  a  28-Foot  Draft 


Figure  24  -  Residual  and  Wavemaking  Resistance  Coefficients  versus  Speed-Length  Ratio 
for  SWATH  IV  (Model  5287)  at  a  25-Foot  Draft 
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Figure  33  -  Residual  and  Wavemaking  Resistance  Coefficients  versus  Speed-Length  Ratio 
for  SWATH  V  (Model  5301-1)  at  a  32-Foot  Draft 


MODEL 
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Figure  36  -  Residual  and  Wavemaking  Resistance  Coefficients  versus  Speed-Length  Ratio 
for  SWATH  &C-2  at  a  28-Foot  Draft 
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